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ABSTRACT 

A  six  degree  of  freedom  mathematical  model  of  an  assault 
boat  landing  on  and  retracting  off  a  sloping  beach  through  surf 
is  developed  based  on  both  theoretical  and  experimental  results. 
Included  are  the  effects  of  wave  refraction,  winds,  shoaling, 
breaking  and  wave  run-up  on  the  beach.  Control  forces  and 
moments  resulting  from  the  coxswain  ability  to  alter  the  throttles 
and  rudders  are  also  included.  This  report  is  to  serve  as  a 
basis  for  the  development  of  a  real  time  assault  boat  simulator 
for  the  training  of  coxswains. 
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The  present  study  is  to  provide  a  mathematical  representation  of  waves  and 
surf  typical  of  conditions  encountered  in  the  conduct  of  assault  landings. 
The  derived  representations  consist  of  equations  giving  the  varying 
characteristics  of  a  wave  environment  which  affect  the  six  degrees-of- 
freedom  motion  of  an  LCM(6)  vehicle.  These  formulations,  data  and  ideas 
will  be  useful  in  developing  a  real-time  computer  simulation  of  landing 
craft  behavior  in  waves.  The  purpose  of  such  a  simulation  relates  to  the 
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Development  of  the  simulation  at  the  Naval  Training  Device  Center  will 
involve  verification  and  possible  extension  of  the  results  documented  here. 


HARVEYaJ.  saltzman 
Projeoft  Engineer 
Naval  Training  Device  Center 
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Water  depth  at  start  of  sloping  beach 
Moment  of  inertia  about  x,  y,  or  z-axis 
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Wave  amplitude  factor  due  to  refraction 

Moment  about  x.  y.  or  z-axis;  also 
referred  to  as  rolling,  pitching  or 
yawing  moment 

Wave  number 

Wave  number  in  deep  water 

Length  along  load  waterline,  used  as 
reference  length 

Half-length  of  equivalent  rectangular 
block  ship  representation 

Mass 

Local  heave  damping  coefficient 
Propeller  rotational  speed,  rev. /sec. 
Pressure 

Angular  velocity  component  about  x,  y, 
or  z-axis 

Run-up  magnitude 

Spectral  energy  of  waves  referred  to  a 
fixed  position 

Wave  period,  also  used  as  symbol  for 
propeller  thrust 

Time 
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Resultant  wind  speed  relative  to  the 
boat 
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SECTION  I 
INTRODUCTION 

The  development  of  adequate  techniques  for  military  person¬ 
nel  in  carrying  out  their  assigned  tasks  in  connection  with 
vehicle  performance  is  often  arrived  at  by  the  use  of  various 
training  procedures.  For  one  reason  or  another,  it  is  often 
not  possible  or  permissible  to  allow  training  to  be  carried  out 
with  the  use  of  a  vehicle,  whether  it  be  a  ship,  or  aircraft, 
etc.  In  order  to  acquaint  operational  personnel  with  the  ele¬ 
ments  they  will  experience  in  actual  operation  and/or  control  of 
such  vehicles,  recourse  is  made  to  the  use  of  training  devices. 
These  devices  provide  a  simulation  of  the  realistic  environment, 
obtained  as  a  combination  of  visual  displays  and/or  motion  cues. 
In  the  case  of  training  of  an  assault  boat  coxswain,  there 
presently  exists  a  limitation  due  to  the  unavailability  of 
training  boats;  requirements  of  travel  to  training  beach  loca¬ 
tions;  and  possible  adverse  weather  and  surf  conditions.  Thus  a 
need  arises  for  development  of  an  assault  boat  coxswain  trainer, 
where  a  mathematical  model  has  to  be  derived  for  use  in  the 
training  device  computer  in  order  to  represent  the  vehicle  motion 
and  the  wave  conditions. 

This  project  was  carried  out  for  the  Naval  Training  Device 
Center  under  support  of  Contract  No.  N61339-68-C-0152 . 
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SECTION  II 

STATEMENT  OF  THE  PROBLEM 

In  the  development  of  the  trainer,  it  is  necessary  to  repre- 
sent  the  time  history  of  the  waves  that  are  approaching  a  sloping 
beach,  where  the  wave  representation  provides  information  on  the 
wave  system  form  in  the  vicinity  of  the  assault  boa\.  as  it 
\jproaches  the  beach.  The  equations  for  the  waves  will  represent 
the  wave  properties  and  form  up  to  surf  formation  when  breakers 
are  formed  on  the  beach,  and  this  wave  system  will  continually 
affect  the  assault  boat  by  providing  wave-induced  excitations 
that  cause  variations  in  the  six  degrees  of  freedom  of  a  repre¬ 
sentative  assault  boat  such  as  the  LCM  (6) .  In  addition  the 
vehicle  motions  will  be  affected  as  a  result  of  the  coxswain 
control,  which  is  manifested  by  his  ability  to  control  the 
throttle,  gear  shift,  and  rudder,  tnereby  determining  the  forward 
speed  and  the  ship  heading.  The  effects  of  wind  are  also  included, 
since  strong  winds  will  create  forces  on  the  hull  and  will  affect 
its  motion  in  the  horizontal  plane,  i.e.  in  the  axial  and  lateral 
directions . 

The  present  report  presents  the  results  of  a  combined  theo¬ 
retical  and  experimental  study  which  has  as  its  objective  the 
development  of  a  mathematical  model  that  represents  the  six 
degrees  of  freedom  motion  of  assault  boats  in  waves.  The  model 
will  a.i-o  provide  a  time  history  of  the  waves  along  a  sloping 
beach,  in  proximity  to  the  assault  boat,  and  it  will  allow  treat¬ 
ment  of  the  approach  to  the  beach  as  well  as  the  retraction 
operation  of  the  boat  returning  from  the  shore  to  seaward.  An 
important  consideration  in  the  mathematical  model  is  the  ability 
to  obtain  real  time  computer  solutions,  thereby  permitting 
direct  utilization  in  a  training  device  with  a  man  "in  the  loop" 
to  provide  control  actions  and  observe  their  consequences.  The 
mathematical  model  is  presented  in  the  present  report,  based 
upon  many  results  obtained  in  an  experimental  model  study  in  a 
towing  tank  [l].The  data  analyses,  approximations,  computations, 
etc.  providing  further  support  for  the  proposed  mathematical 
model,  including  verification  and  validation  studies,  are 
presented  in  an  associated  report  [2]  that  represents  the  final 
descriptive  report  on  this  project. 
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SECTION  III 
METHOD  OF  PROCEDURE 


A.  ANALYSIS 

1.  GENERAL 

The  techniques  used  in  deriving  the  equations  presented 
herein  represent  the  application  of  the  present  state  of  the  art 
in  accounting  for  motions  of  ships  in  waves.  As  such  linear 
theory  is  primarily  used,  with  appropriate  expansion  to  encompass 
certain  nonlinear  effects  when  they  are  considered  to  be  signifi¬ 
cant.  A  similar  treatment  is  applied  when  considering  the  waves 
as  they  progress  into  shore  and  affect  the  boat  motion.  However 
since  the  major  objective  of  this  study  is  to  obtain  a  mathemati¬ 
cal  model  that  will  suffice  for  application  in  a  training  device, 
certain  approximations  are  made  in  order  to  allow  simplification 
and  achieve  results  useful  for  that  purpose.  Thus  the  precision 
associated  with  a  detailed  research  program  directed  toward 
shallow  water  wave  effects  on  ships;  complete  details  of  wave 
form  for  breakers;  determination  of  the  change  in  all  of  the 
spectral  properties  of  waves  along  a  sloping  beach;  etc.  is  not 
achieved  in  this  work,  but  nevertheless  a  utilitarian  tool 
adequate  for  training  purposes  is  obtained.  The  foregoing  repre¬ 
sents  the  fundamental  philosophy  underlying  the  development  of 
the  mathematical  model  described  in  this  report. 

In  formulating  the  equations  of  motion  all  six  degrees 
of  freedom  of  the  boat  are  represented.  In  the  usual  case  of  a 
ship  in  waves,  the  vertical  plane  motions  of  pitch  and  heave  are 
formulated  according  to  linear  theory  (see  [3J  and  [4])  for  a  fixed 
forward  speed,  and  the  effect  of  surge  is  neglected.  However  in 
the  present  case  where  forward  speed  control  is  manifested  by 
the  coxswain  controlling  the  engine,  and  also  using  the  ability 
to  apply  differential  thrust  as  a  means  of  controlling  heading, 
the  surge  degree  of  freedom  (i.e.  axial  motion)  is  coupled  with 
the  lateral  plane  motions  as  well  as  the  vertical  plane  motions. 

The  nature  of  the  coupling  is  both  linear  and  nonlinear,  and  the 
functional  form  resulting  will  be  the  most  simple  possible  in 
order  to  achieve  the  required  ease  of  computer  representation 
and  soluv.ion. 

•Jhe  wave  motion  is  represented  in  a  time  history  form 
for  regions  in  proximity  to  the  boat,  thereby  allowing  incorpor¬ 
ation  into  a  proposed  training  device  since  the  view  of  the  local 
wave  environment  is  then  available  in  conjunction  with  a  simu¬ 
lated  boat  motion.  The  wave  representation  includes  the  effects 
of  shoaling,  which  represents  the  influence  of  the  sloping  beach 
depth  changes  on  the  wave  properties,  as  well  as  refraction 
effects  that  account  for  the  influence  of  changing  depth  on  the 
direction  of  wave  propagation  as  the  shoreline  is  approached. 
Methods  allowing  evaluation  of  deterministic  properties  of 
waves  of  sinusoidal  form  (initially,  in  deep  water)  as  they 
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progress  onto  the  beach,  as  well  as  for  a  prescribed  spectral 
distribution  of  irregular  random  waves,  are  included. 

The  equations  of  motion  for  the  boat  in  regular  waves 
are  formulated  on  the  basis  of  linear  theory  initially,  by  the 
balance  of  inertial,  damping,  restoring,  coupling,  and  wave 
excitation  forces  and  moments.  Both  hydrodynamic  and  hydro- 
static  effects  of  body-fluid  interaction  are  included,  with 
many  terms  determined  by  application  of  the  technique  denoted 
as  strip  theory  (a  form  of  s lender-body  theory) ,  as  demonstrated 
in  [3]  and  [4] .  This  method  makes  the  assumption  that,  for  an 
elongated  body  with  transverse  dimensions  small  compared  to  its 
length,  the  flow  at  any  cross-section  is  independent  of  the  flow 
at  any  other  section,  and  hence  the  flow  problem  is  reduced  to  a 
two-dimensional  problem  in  the  transverse  plane.  The  forces  at 
each  section  are  found  by  this  method,  and  the  total  force  is 
found  by  integrating  over  the  body  length.  While  this  theo¬ 
retical  technique  has  been  found  useful  for  heave  and  pitch 
motions  in  waves  for  a  ship  at  a  fixed  forward  speed,  the 
presence  of  the  other  m^aes  of  motion  together  with  the  influ¬ 
ence  of  forward  speed  control  requires  obtaining  extensive  model 
test  data  for  the  various  hydrodynamic  forces  that  act  on  a 
landing  cr  ft.  Accordingly  a  model  test  program  was  carried  out 
for  this  purpose,  and  the  results  presented  in  [1].  The  data  from 
that  report,  when  analyzed  and  reduced  to  the  required  final  form, 
yielded  much  of  the  information  required  for  application  to  the 
present  mathematical  model,  and  a  description  of  the  techniques 
used  to  obtain  the  necessary  hydrodynamic  information  is  presented 
in  [2] ,  the  complete  final  report  for  this  project. 

The  particular  assault  boat  studied  in  detail  in  this 
program  is  the  LCM  (6).  Information  on  the  operation  of  this 
craft,  as  well  as  some  description  of  its  physical  character¬ 
istics,  is  provided  in  [5],  Body  lines,  section  characteristics, 
rudder  profile,  propeller  outlines,  etc.  were  obtained  from  [6]. 
This  particular  craft  has  twin  screws,  both  of  which  are  right- 
handed  and  offset  from  the  ship  centerline.  Twin  rudders  are 
also  present,  protected  by  skegs,  and  lying  aft  of  the  propellers. 
The  rudders  operate  in  unison,  and  they  are  controlled  by  one 
wheel  in  the  pilothouse.  Two  displacements  are  considered,  the 
heavier  one  being  the  condition  when  the  craft  approaches  the 
beach  fully  loaded  and  the  lighter  being  the  condition  for  the 
craft  leaving  the  beach  empty.  A  summary  of  the  numerical  values 
of  the  important  parameters  characterizing  the  LCM  (6) ,  for  which 
the  present  mathematical  model  was  derived,  is  presented  in 
Table  1. 


2.  WAVE  CHARACTERISTICS  IN  SHALLOW  W?  T*ER 

In  order  to  illustrate  the  basic  characteristics  of 
waves  as  they  propagate  toward  a  beach  from  deep  water,  the 
fundamental  equations  for  amplitude,  wavelanatn,  and  direction 
changes  are  presented  in  the  follow  ny  T.iitially  there  are 
considered  to  be  plane  waves  offshore  o<'  frequency  u  (rad./sec.), 
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amplitude  aQ,  and  wave  direction  angle  and  they  are  travelling 

shoreward  over  a  prismatic  bottom  towards  a  parallel  shoreline, 
as  shown  in  Figure  1.  With  the  assumption  of  conservation  of 
energy  and  frequency  and  no  reflection  of  the  waves,  as  in  the 
case  of  a  slowly  varying  bottom  contour  of  known  depth  h(x), 
the  inshore  wavelength  A,  amplitude  a  and  wave  direction  angle 
e  can  be  predicted;  however  no  attempt  is  made  in  the  analysis 
to  determine  the  phase  of  the  waves.  Assuming  small  amplitude 
and  using  the  principle  of  linear  superposition  this  can  also 
be  extended  to  an  offshore  spectrum  of  waves  with  amplitudes 
distributed  according  to  frequency. 


For  the  case  of  a  simple  plane  wave  of  frequency  to 
travelling  in  a  direction  perpendicular  to  the  beach,  hydro 
dynamic  wave  theory  £71 yields  the  result 


where 


and 


a  =  aQ  D (to) 


D  = 


1  + 


2kh 

sTnH'"2Jch' 


k_  =  =  k  tanh  kh 


(1) 

(2) 

(3) 


A  graph  of  the  relationship  in  Equation  (3)  is  shown  in  Figure  2, 
together  with  other  relations  covering  wavelength,  slope,  etc. 

In  all  cases  inshore  k  >  k  is  valid,  so  that  the  corresponding 

wavelengths  A  *  ate  therefore  less  than  those  offshore.  The 

wave  amplitudes  are  found  to  first  decrease  in  moderately  deep 
water  and  then  increase  in  very  shallow  water,  but  in  all  cases 
the  wave  steepness  increases  in  shallow  water.  All  of  the  results 
obtained  in  the  foregoing  analysis  are  general,  and  not  restricted 
to  any  particular  bottom  contour. 

In  the  case  of  waves  approaching  the  beach  at  an  offshore 
angle  to  the  beach  6Q,  refraction  takes  place.  For  a  beach  whose 

depth  decreases  monotonically  with  x  inshore,  the  wave  angles 
decrease  and  the  crests  tend  to  be  parallel  to  the  beach  far 
inshore.  The  controlling  relationship,  referred  to  as  "Snells 
Law"  [7]  is  k 

sin  6  =  sin  (4) 


There  are  corresponding  amplitude  variations  due  to  converging 
of  the  wave  rays,  so  that  in  the  case  of  refraction 

a  =  aoKD  ,  (5) 
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where 


(6) 


and  D  is  as  given  previously  in  Equation  (2) .  An  illustration 
of  the  refraction  of  waves  is  given  in  Figure  3,  for  the  case 
where  =  45°  and  the  beach  is  straight  with  parallel  depth 

contours.  The  trend  toward  crests  parallel  with  the  beach 
will  maintain  itself  until  breaking  occurs,  which  will  determine 
the  actual  direction  of  the  waves  relative  to  the  beach  at 
the  shoreline. 

The  wave  system  characteristics  described  above  are 
based  on  linear  theory,  and  hence  have  a  restricted  validity 
as  the  wave  progresses  into  more  shallow  water  and  approaches 
the  shoreline.  The  wave  shape  will  begin  to  change  from  the 
sinusoidal  form  and  the  effects  of  breaking  then  manifest 
themselves,  with  the  phenomenon  then  being  nonlinear  and  complex. 
A  discussion  of  these  effects  will  be  given  in  a  later  section, 
while  the  present  discussion  is  sufficient  to  provide  a  general 
description  of  the  major  overall  properties  of  the  waves. 

A  deterministic  approach,  similar  to  that  in  [8]  and 
[9],  can  be  applied  to  the  case  of  waves  travelling  into  shallow 
water  from  an  offshore  deep  depth  location.  In  order  to  simplify 
the  analysis  the  waves  offshore  are  assumed  at  first  to  be 
travelling  perpendicular  to  the  beach,  and  the  waves  are  measured 
with  respect  to  an  absolute  frame  of  reference  rather  than  with 
respect  to  the  moving  craft. 

Consider  waves  as  shown  in  Figure  4  moving  from  offshore 
to  inshore  in  the  direction  of  x  over  a  bottom  of  slowly  varying 
depth  h (x) ,  with  the  latter  variation  slow  enough  to  insure 
complete  transmission  of  energy  without  reflection.  This 
results  in  the  connection  between  the  amplitude  a  of  the  waves 
and  the  offshore  amplitude  a  ,  discussed  previously  and  defined 
by  Equations  (l)-{3). 

A  precise  presentation  of  the  waves  at  a  particular 
point,  with  reference  to  their  form  (as  a  function  of  time)  at 
another  point,  can  be  achieved  in  a  simple  manner,  where  any 
change  in  phase  is  indicated  directly.  This  can  be  done  for 
a  single  sinusoidal  wave  form  and  generalized  to  an  arbitrary 
wave  system,  following  the  techniques  in  [8],  [9].  The  same 
type  of  result  can  also  be  obtained  for  the  case  where  the 
waves  travel  to  some  angle  other  than  normal  to  the  beach, 
which  would  include  the  effect  of  refraction  as  well.  However 
a  detailed  consideration  of  this  particular  technique  was  found 
to  be  difficult  to  implement  for  the  present  problem  due  to  the 
requirement  of  accounting  for  the  influence  of  the  craft  forward 
speed  (which  varies  with  time) ;  anomalies  associated  with 
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relations  between  wave  frequencies  and  frequency  of  encounter  in 
following  seas;  and  expected  computational  requirements  that 
inhibit  real-time  simulation  of  the  wave  system  in  this  case. 

This  approach  would  only  have  use  for  determining  the  exact 
representation  of  a  simplified  wave  system ,  or  for  a  spectral 
distribution  of  irregular  waves,  if  the  instantaneous  phase 
were  required.  Since  the  application  of  the  results  in  this 
study  are  for  simulation  purposes,  only  conditions  representative 
of  what  occurs  are  required.  Actual  wave  systems  at  sea  are 
irregular  random  processes  [10],  and  they  are  represented  as  a 
summation  of  many  individual  waves  of  various  frequencies,  each 
with  a  random  phase,  so  that  the  necessity  for  representing 
phase  relationships  precisely  is  not  present  for  this  investiga¬ 
tion.  Thus  a  simplified  treatment  for  representing  component 
waves  that  make  up  an  actual  wave  spectrum  is  used,  thereby 
allowing  ease  of  computation  as  well  as  a  simple  means  of  repre¬ 
senting  the  wave  forces  associated  with  such  wave  systems. 

To  illustrate  the  above,  consider  a  single  component 
plane  wave  of  frequency  to  and  amplitude  aQ  offshore  propagating 

inshore  towards  a  beach  located  at  x  =  x^.  The  special  case 

of  a  constant  slope  linear  beach  defined  by 

h  =  h  -  ax 
o 

with 


is  the  prescribed  depth  variation  considered  throughout  this 
study,  for  simplicity.  The  offshore  wave  is  represented  by 

n  =  aQ  sin  (kx-ust)  ,  (9) 

and  at  the  point  x^  the  wave  is  given  by 

n  =  a-^  sin  (k^x^-wt)  (10) 

where 


(7) 

(8) 


al  ■  ao  D1 


with  the  values  of  and  k^  found  from  Equations  (2)  and  (3) , 
using  the  local  depth  value 


*  hQ  -  a  x1  (12) 

Similarly  at  the  point  x2 , 

n  -  a„  sin  (k0x--wt)  (13) 

i.  4 

with  similar  relations  for  a2,  k2/  and  D2  in  terms  of  the  local 
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depth  etc. 

In  order  to  represent  the  effect  of  refraction  on  the 
preceding  results,  a  single  component  plane  wave  of  frequency 
to  at  an  initial  angle  Bc  offshore  is  considered  progressing 

inshore.  Assuming  no  reflection  from  the  beach,  the  wave 
continuously  refracts  at  the  angle  B  until  it  finally  arrives 
at  the  beach,  where  the  relationship  between  B  and  BQ  is 

given  by  Snell's  law,  Equation  (4).  The  amplitude  of  a  plane 
wave  travelling  inshore  at  an  angle  B  to  the  x-axis  is 


=  a  sin  k(x  cos  B  +  y  sin  8)  ~  wt 

L.  <Jl 

=  a  sin . |kxx  +  k^y  -  wtj 

/ 

(14) 

where 

k  =  k  cos  B 

X  (15) 

k  =  k  sin  B  =  k  sin  B^  =  constant 
y  o  o 

(according  to  Snell's  law). 

Refraction  and  shoaling  affect  the  wave  amplitude 
according  to 


a.  =  K.D.art 
1  1  1  o 


(16) 


where  and  are  defined  in  Equations  (2)  and  (6),  and  the 

i-subscript  relates  the  spatial  locations  considered.  The 
latter  relation  assumes  no  reflection  of  energy  while  proceeding 
into  the  beach  and  total  absorption  of  energy  at  the  beach,  as 
noted  previously.  The  relations  above  provide  the  means  to 
predict  the  detailed  behavior  of  waves  at  different  inshore 
points  given  the  beacn  slope  and  the  frequency,  amplitude  and 
direction  of  the  waves  at  an  offshore  point.  Thus, 


nl 


a^  cos 


(17) 


at  point  (xi»y^)»  and 


02  ~  &2  cos 


x0  +  k 

2  y- 


(18) 


at  (x2,y2) •  in  accordance  with  the  method  described  above,  where 


ai 


K2D2 

K1D1 


(19) 
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and  the  K. -values  are  found  in  terms  of  the  local  depth  h.. 

The  preceding  discussion  considered  the  characteristics 
of  the  waves  when  viewed  from  an  absolute  frame  of  reference,  but 
the  present  problem  requires  consideration  of  the  wave  character¬ 
istics,  time  history,  etc.  in  proximity  to  the  assault  boat  while 
it  is  moving  at  some  forward  speed.  When  the  boat  approaches  a 
shoreline,  the  seas  are  generally  in  the  same  direction  as  the 
craft  forward  motion,  and  hence  are  following  and  overtaking  seas. 
When  the  craft  is  retracting  from  a  beach  and  going  seaward  it 
will  be  proceeding  opposite  to  the  wave  propagation  direction, 
and  hence  the  seas  will  be  head  seas  relative  to  the  boat.  The 
main  effect  of  the  forward  speed  of  the  boat  is  to  alter  the 
frequency  of  encounter  of  the  waves  relative  to  the  boat,  with 
different  influence  and  consequences  in  each  case.  Assuming  pure 
head  sea  operation,  the  frequency  of  encounter  in  deep  water  is 
given  by  2 

10  =  to  4-  u  ,  (20) 

®  g 


assuming  constant  forward  speed  u.  For  the  situation  where  the 
boat  approaches  the  beach,  with  the  waves  as  pure  following  seas, 
the  frequency  of  encounter  (assuming  deep  water  and  constant 
speed)  is  represented  by 

..2., 

(21) 


=  to  - 


0)2U 


The  change  in  the  frequency  of  encounter  due  to  forward 
speed  in  shallow  water,  just  as  in  the  case  of  deep  water,  is 
exhibited  by  the  following  relations: 

w  «  U)  +  ku  (22) 

e 

for  head  seas,  and 

Hi.  =  iii  -  ku  (23) 

e 

for  following  seas.  Since  the  wave  number  k  varies  with  position 
along  the  beach  in  the  shallow  water  case,  the  frequency  of  encounter 
will  vary  as  the  boat  moves  along  the  sloping  beach.  However,  a 
knowledge  of  the  ship  speed  and  its  position  along  the  beach 
(thereby  providing  information  on  the  local  depth)  will  allow 
computation  of  the  change  in  frequency  of  encounter. 

The  effect  of  forward  speed  also  manifests  itself  ir.  ?. 
reshaping  of  the  wave  spectrum  as  a  function  of  (see  (10) ). 

This  effect  is  important  if  a  continuous  change  in  the  wave 
spectrum  is  to  be  applied  as  a  function  of  craft  position  and 
time,  which  would  be  an  impractical  technique  for  simulation 
purposes.  Consideration  of  this  encounter  spectrum  only  provides 
insight  into  the  particular  wavelengths  that  would  have  signifi¬ 
cant  influence  on  ship  motions,  but  J  sat  is  only  important  f >r 
purposes  of  analysis  rather  than  application  to  simulation,  which 
is  the  prime  purpose  of  this  study. 
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Since  the  forward  speed  of  the  craft  will  generally 
vary  due  to  wave  effects  and  coxswain  control  inputs,  the  use  of 
a  frequency  of  encounter  is  only  limited  to  interpretive  purposes 
and  simplified  illustrations.  When  considering  wave  representations, 
wave-induced  forces,  etc,  relative  to  the  moving  craft  the  "argu¬ 
ment"  of  the  defining  sinusoidal  functions  will  be 


•X 

+1 

+1 

3 

u  dt 

\  J 

and  discussion  of  the  applicability  of  this  term  will  be  provided 
in  a  later  section. 

3.  BREAKER  AND  SURF  CHARACTERISTICS 

As  indicated  previously,  the  wave,  height  of  a  wave 
approaching  a  shoreline  will  increase  as  it  progresses  closer 
to  the  shore,  and  eventually  reaches  a  condition  where  the  wave 
form  changes  and  the  wave  then  becomes  unstable  and  breaks.  The 
breaking  is  characterized  by  the  formation  of  turbulence  as  water 
spills  from  the  crest  while  the  wave  advances,  and  also  by  means 
of  the  entire  crest  falling  forward  in  a  plunging  form  and 
dissipating  its  energy  in  the  turbulence  associated  with  the 
breaking  process.  The  surf  zone  is  a  region  where  the  waves 
tend  to  be  quite  steep,  and  it  extends  from  the  outermost 
breakers  to  the  region  where  the  water  is  rushing  up  the  beach. 

In  this  region  nonlinear  processes  become  very  important,  and  con¬ 
sequently  it  is  the  region  which  is  least  amenable  to  theoretical 
treatment.  Many  observations  of  breaker  characteristics  have 
been  made  experimentally,  and  the  nature  of  the  breakers  is 
defined  in  accordance  with  the  manner  of  breaking,  as  a  function 
of  the  original  deep  water  wave  slope  and  the  beach  slope  (see 
[11]). 

When  regular  waves  propagate  inshore,  they  retain  their 
sinusoidal  outline  until  the  ratio  of  the  depth  to  deep  water 
wavelength  reaches  the  condition 

h  Hn  H 

~  =  .04  +  1,35  ~  ,  :or  t-2  >  .01  (25) 

o  o  K  o 

where  Hq  and  \Q  are  the  deep  water  height  and  wavelength.  This 

relation  is  based  upon  results  presented  in  [11]  ,  where  the 
minimum  depth  to  which  wave  height  theories  apply  for  wave 
transformations  in  shoaling  water  is  illustrated.  After  that 
point  considerable  departure  from  the  theoretical  outline  occurs, 
with  such  departures  dependent  on  the  original  wave  steepness. 

The  hydrodynamic  theory  of  waves  of  finite  height  [7]  indicates 
that  breaking  occurs  when  the  crest  angle  is  120°,  with  the 
particle  speed  at  the  crest  equal  to  the  wave  speed  c,  and  the 
particle  acceleration  equal  to  tl  at  of  gravity.  The  maximum 
steepness  for  such  a  wave  in  deep  water  is  found  to  be 
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H  1 

-  =  0.142  »  y  ,  with  an  actual  indication  of  the  phase  speed 

about  10%  greater  than  that  from  linear  theory.  Since  the 
breaking  of  waves  is  a  nonlinear  feature,  only  limited  informa¬ 
tion  for  mathematical  models  can  be  obtained  from  a  linear 
analysis,  and  an  approximation  is  necessary  in  the  present  study. 


The  wave  representation  in  the  present  investigation 
will  not  allow  "foaming",  "white  horses",  etc.  that  dissipate 
wave  energy  by  partial  breaking  at  the  top  of  the  crest.  This 
is  usually  associated  with  steep  waves  in  deeper  water  as  they 
run  into  shallow  water,  and  they  break  so  that  the  maximum  steep¬ 
ness  is  maintained  at  a  certain  level.  The  breaking  character¬ 
istics  for  such  waves  are  quite  different  from  those  for  low 
long  waves  offshore,  which  are  a  more  appropriate  condition  for 
the  present  case.  For  those  waves  -the  wave  profile  inshore 
tends  to  concentrate  the  wave  energy  into  the  crest,  where  the 
crest  rises  and  occupies  only  a  small  fraction  of  the  total 
wavelength,  with  the  trough  becoming  relatively  flat  and 
occupying  the  major  portion  of  the  wavelength.  Thus  the  final 
pattern  is  a  series  of  isolated  steep  crests  separated  from 
each  other  by  long  flat  troughs  which  appear  to  have  only  a 
little  depression  below  the  mean  water  level.  Such  waves  then 
break  at  relatively  large  heights  producing  a  significant 
"run-up"  along  the  beach. 

An  approximation  to  the  wave  characteristics  in  shallow 
water  just  prior  to  breaking  cam  be  obtained  by  use  of  the  soli¬ 
tary  wave  concept.  While  the  agreement  is  not  precise  for 
certain  aspects  of  the  wave  characteristics,  the  approximation 
of  the  solitary  wave  is  a  sufficiently  good  representation  that 
allows  a  useful  mathematical  form  for  the  waves.  This  form  is 
also  applicable  to  the  determination  of  the  excitation  forces 
associated  with  such  a  wave  acting  on  a  boat  as  it  approaches 
the  beach.  Thus  a  mathematical  model  of  the  solitary  wave  is 
used  for  such  purposes  in  this  study. 

Since  the  breaker  height  and  the  depth  at  which  breaking 
occurs  depend  upon  the  beach  slope  and  the  deep  water  wave  slope, 
a  single  relationship  that  is  unique  for  all  conditions  is 
difficult  to  find.  However,  as  an  approximation,  the  limiting 
condition  for  solitary  wave  height  relative  to  the  depth  is  used 
as  a  representation  of  an  effective  criterion  for  establishment 
of  the  breaker  conditions.  This  result  is  represented  by  the 
relation 


%  5 


«  .714 


(26) 


where  the  b-subscript  corresponds  to  breaking  conditions,  which 
will  be  used  as  an  approximation  in  the  present  work.  The 
procedure  to  be  used  in  the  mathematical  model  is  to  determine 
the  local  (maximum)  wave  height  at  the  location  corresponding 
to  the  condition  given  in  Equation  (25),  and  then  to  assume 
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that  the  solitary  wave  has  the  height  corresponding  to  that 
magnitude  at  that  point.  The  solitary  wave  then  propagates 
along  the  sloping  beach  in  accordance  with  the  basic  relations 
defining  such  a  wave.  A  series  of  solitary  waves  can  repeat 
themselves  at  the  same  period  as  a  single  oncoming  sinusoidal- 
form  wave  assumed  present  in  deep  water  offshore  (the  original 
wave) ,  but  in  the  case  of  a  spectral  distribution  of  waves  off¬ 
shore  an  approximation  can  be  made  to  use  the  period  of  maximum 
spectral  energy  of  the  particular  deep  water  spectral  distribu¬ 
tion. 


The  height  of  each  particular  solitary  wave  will  vary 
with  the  local  depth  in  the  same  manner  as  an  ordinary  small 
amplitude  wave,  i.e.  according  to  the  classical  Green's  law 
for  shoaling  waves,  i.e. 


(27) 


which  relates  the  heights  at  two  points  in  terms  of  their  local 
depths.  The  solitary  wave  form  propagates  in  this  manner  until 
the  limit  condition  given  by  Equation  (26)  is  achieved,  at  which 
time  breaking  can  be  said  to  occur.  The  basis  for  using  the 
solitary  wave  is  the  similar  appearance  of  solitary  waves  to 
breakers  (i.e.  spilling  breakers)  and  the  tendency  to  form  a 
relatively  high  /pulse-like  wave  elevation.  The  approach  is  that 
of  a  quasi-steady  variation,  where  at  each  instant  a  steady 
state  representation  of  the  solitary  wave  is  used  in  accordance 
with  the  variations  caused  by  local  depth  changes  along  the 
sloping  beach.  A  practical  reason  for  use  of  this  model  is  the 
complicated  requirement  of  applying  numerical  techniques  such 
as  the  method  of  characteristics  (e.g.  [12])  or  the  Lagrangian 
method  of  [13]  to  continuously  compute  the  wave  form  as  an  actual 
unsteady  problem  up  to  the  point  of  breaking.  The  actual  breaking 
conditions  are  not  properly  represented  by  either  of  the  latter 
approaches,  and  their  computational  complexity  limits  their 
utility.  This  is  especially  significant  since  it  is  necessary 
x  *'  use  the  hydrodynamic  wave  system  to  determine  the  wave-induced 
forces  on  a  craft  in  a  real-time  computation  where  representation 
of  the  waves  in  both  space  and  time  is  required. 


The  equation  representing  the  wave  form  for  the  solitary 
wave,  assuming  that  it  is  propagating  along  the  positive  x^- 
direction,  i3  given  by  [11]  as  _ 


=  H  sech2 


where 


yir  (*0-c« 


(28) 


c  =  /5h  |i  +  |  2 


*  /g (h+H) 


(291 


is  the  wave  propagation  speed,  and  H  is  the  height  of  the  wave 
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above  the  still-water  level  corresponding  to  the  local  depth  h. 
The  expression  in  Equation  (28)  provides  a  representation  of  the 
variation  of  the  wave  form  as  a  function  of  space  and  time,  where 
the  height  changes  in  accordance  with  the  relation  given  in 
Equation  (27)  as  the  wave  progresses  along  the  beach.  This 
expression  is  valid  for  the  case  of  the  wave  propagating  in  a 
direction  normal  to  the  shoreline,  as  observed  in  a  fixed  frame 
of  reference. 


Since  use  of  the  solitary  wave  representation  as  an 
equivalent  to  the  breaker  form  oc  curs  in  relatively  shallow 
depths,  it  is  assumed  that  no  further  refraction  occurs  and  that 
the  wave  propagates  at  the  angle  appropriate  to  the  condition 
where  "transition"  occurs  between  the  sinusoidal  wave  forms  and 
the  solitary  wave.  The  representation  of  the  wave,  at  the  angle 
corresponding  to  the  transition  orientation  (denoted  as  3fc) ,  is 
given  by 


nm  =  H  sech2 


L 


3H_ 
4h3 


x  cos  6.  +  y^sxn  8,  -  ct 
o  t  o  t 


jj 


(30) 


This  wave  orientation  will  maintain  itself  up  to  breaking,  beyond 
which  no  further  wave  form  is  considered  to  exist. 


Following  the  occurrence  of  breaking,  a  form  of  "wave  of 
translation"  or  bore  is  assumed  to  exist  in  the  water,  and  this 
produces  a  certain  run-up  on  the  beach,  where  the  run-up  is  the 
vertical  distance  between  the  maximum  height  the  water  runs  up 
the  beach  and  the  still-water  level.  The  magnitude  of  the  run-up 
is  given  by  [14],  as  obtained  from  a  large  series  of  experiments. 

It  is  given  as  a  function  of  the  deep  water  wave  height  and  period, 
and  for  the  present  investigation  the  deep  water  wave  height 
reference  is  selected  as  the  significant  height  offshore,  where 
the  significant  height  is  the  mean  of  the  one-third  highest  waves, 
which  is  given  by 


H 


o 


sig. 


(31) 


with  o  the  rms  wave  amplitude  in  deep  water.  The  reference 


period  for  the  evaluation  of  the  run-up  is  the  period  of  maximum 


2  77 

spectral  energy,  T  =  —  .  In  the  evaluation  of  the  conditions  for 

wm 

the  transition  to  a  solitary  wave,  expressed  by  Equation  (25), 
the  deep  water  wavelength  chosen  for  that  case  is  the  value 
associated  with  the  deep  water  wave  frequency  m  ,  co^c' stent 


with  the  assumptions  used  in  this  treatment. 


It  is  very  difficult  to  predirt  t  magnitude  of  the 
velocities  in  the  water  advancing  up  ti  *  each  following  breaking, 
since  that  is  dependent  on  the  value  of  the  beach  slope,  friction 
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effects  on  the  bottom,  fluid  mixing  processes,  etc.  However,  it 
is  possible  to  approximate  the  value  of  the  longitudinal  veloci¬ 
ties  in  the  fluid  after  breaking  by  applying  some  of  the  results 
of  [15] ,  with  the  assumption  that  the  water  region  after  breaking 
is  represented  by  a  constant  increased  height  about  the  undis¬ 
turbed  water  level  equal  to  the  run-up  R.  Thus  a  uniform  height 
of  fluid  is  assumed  to  exist  from  the  end  of  the  breaking  wave  ip 
to  the  shoreline,  where  the  instantaneous  height  of  water  above 
the  "ground"  level  of  the  beach  is  (h+R) ,  where  h  is  the  local 
depth  below  the  undisturbed  water  level.  Assuming  that  tire 
propagation  speed  of  this  bore  is 

cR  =  /glh+R)  ,  (32) 

similar  to  the  relation  for  the  solitary  wave  given  in  Equation 
(29),  the  results  in  (153  indicate  that  the  longitudinal  particle 
velocities  in  this  bore  are  represented  by 

uR  «  .5  cR  (33) 

This  rush  of  water  toward  shore  is  assumed  to  occur  after  breaking, 

T  2tt 

and  exists  for  a  period  of  about  j  ,  where  T  -  ~  ,  and  this 

m 

repeats  itself  after  each  breaking  wave. 

From  the  point  of  view  of  action  on  a  boat,  when  it  is 
landing,  stationary  on  the  beach,  or  retracting  (as  long  as  it 
is  located  in  the  region  between  the  breaking  wave  and  the 
shoreline) ,  the  effect  of  this  fluid  motion  is  also  to  elevate 
the  fluid  level  by  the  magnitude  of  run-up  R  above  the  still 
water  level.  This  change  occurs  as  a  time-varying  form  which 
can  be  represented  initially  as  a  rectangular  pulse  with  some 
sort  of  decay  so  that  the  return  to  level  water  conditions  is 
present.  Since  it  is  difficult  to  deal  precisely  with  the 
seaward  return  and  backwash  of  the  water  mass  following  breaking, 
this  method  of  approximation  is  considered  sufficient.  Techniques 
for  carrying  out  the  necessary  computations  (with  a  digital 
computer)  that  represent  the  various  aspects  of  the  waves,  viz. 
oscillatory  wave  systems,  solitary  waves,  and  the  bore  character¬ 
istics  following  breaking,  are  discussed  in  the  following 
section  of  the  report. 

4  *  METHODS  FOR  SIMULATION  OF  WAVE  PATTERNS 

In  order  to  represent  the  waves  at  various  points  in¬ 
shore  as  the  boat  progresses,  it  is  necessary  to  determine 
them  from  the  initial  deep  water  waver*  that  exist  further  out 
at  sea.  For  the  present  case  the  waves  must  be  determined  as 
functions  of  time  and  spatial  location,  whether  observed  at  a 
fixed  point  or  when  located  on  a  moving  craft.  A  simulation  tool 
that  has  been  used  extensively  in  representing  waves,  when 
considering  a  random  sea,  is  to  obtain  a  time  history  record 
from  the  output  of  a  white  noise  generator  that  is  sent 
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through  a  spectrum-shaping  filter.  The  amplitude  characteristics 
of  this  filter  are  proportional  to  the  square  root  of  the  desired 
wave  power  spectral  density  at  each  frequency  in  the  wave  band¬ 
width,  and  the  output  signal  is  then  a  single  realization  of  a 
time  series  record  having  the  desired  power  spectrum.  This 
record  serves  as  the  deep  water  wave  system,  and  appropriate 
operation  upon  it  can  produce  the  inshore  waves  and  records  of 
other  wave-related  quantities  via  convolution  methods  [8]  (when 
that  method  is  applicable)  or  by  varying  the  spectrum  in  terms 
of  local  depth.  An  alternate  method  can  be  applied  to  a  given 
offshore  wave  spectrum,  where  the  spectrum  is  divided  into  a 
number  of  strips  of  different  bandwidths  and  each  strip  is  repre¬ 
sented  by  a  sinusoid  whose  magnitude  and  frequency  are  given  by 
the  square  root  of  the  area  of  the  strip  and  by  the  median 
frequency  of  the  strip.  This  representation  produces  a  wave 
output  record  that  appears  random  and  is  a  sufficient  basis  for 
representation  of  the  particular  wave  spectrum  of  interest.  When 
the  spectrum  is  divided  into  a  finite  number  of  sinusoids,  it  is 
relatively  easy  to  determine  the  resultant  inshore  waves,  wave 
forces,  and  other  quantities  used  in  a  simulation  study.  The 
number  of  waves  chosen  to  represent  a  given  offshore  spectrum, 
and  an  illustration  for  the  particular  case  of  Sea  State  4, 
using  8  component  waves,  is  shown  in  Figure  5. 

Each  of  the  constituent  sinusoidal  functions  in  deep 
water  represents  the  far  offshore  waves,  and  they  are  operated 
upon  in  a  computer  by  factors  that  account  for  the  effects  of 
changing  depth  and  rei- action  as  the  wave  proceeds  along  a 
prescribed  sloping  beach,  where  the  depth  is  known  as  a  function 
of  distance  along  the  beach.  The  sum  of  all  the  constituent 
wave  elements,  each  of  which  is  treated  by  the  parameters  changing 
the  amplitude,  wave  number,  etc.  as  indicated  previously,  will 
then  yield  the  complete  wave  form  at  various  points  along  the 
beach.  This  can  also  be  accomplished  for  exhibiting  a  spatial 
variation  at  a  fixed  time,  for  a  series  of  time  increments, 
thereby  providing  the  wave  pattern  in  a  region  as  a  function 
of  time.  This  would  be  the  "picture"  of  the  waves  as  experienced 
by  an  observer  moving  through  the  waves  on  a  boat. 

The  procedure  to  be  used  is  indicated  for  the  case  of  a 
wave  of  frequency  ok  ,  assuming  that  the  wave  propagates  in  a 

direction  normal  to  the  shoreline.  The  wave  is  observed  from  a 
boat  assumed  to  travel  normal  to  the  shoreline,  at  a  speed  u, 
thereby  representing  a  fo  -lowing  sea.  The  origin  of  coordinates 
of  the  inertial  frame  is  located  at  a  point  distant  from  the 

shoreline,  assuming  a  beach  of  constant  slope  a,  and  a  water  depth 
of  hQ  at  this  offshore  reference  point.  With  the  bottom  contour 

given  by 

h  =  hQ  -  a  xQ  ,  (34) 

where  hQ  =  ax^  ,  and  an  initial  wave  elevation  at  that  point 
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(assumed  to  be  in  sufficiently  deep  water)  given  by 

nQ  =  aQ  sin  (k^-o^t)  (35) 

(observed  at  a  stationary  reference) ,  the  wave  as  seen  by  the 
moving  observer  on  the  boat  is  then 


u  dt  -  uj^t  .  (36) 

In  Equation  (36)  the  observation  point  is  chosen  at  a  value  of 
x  measured  relative  to  the  boat,  i.e.  at  a  fixed  distance  from 
the  CG  of  the  craft;  and  other  parameters  such  as  k^  ,  D  are 

dependent  upon  the  local  value  of  h  at  each  instant  of  time,  i.e. 


(37) 


as  the  boat  progresses  toward  the  beach.  The  quantity  k^  is 
defined  by 

“i2 

-  =  k.  tanh  k.h  (38) 

g  l  l 


(as  in  Equation  (3)),  and  0  is  given  by  Equation  (2),  so  that 
the  variation  of  these  parameters  as  the  boat  moves  and  changes 
the  local  value  of  h  is  evident. 


The  process  described  above  is  carried  out  for  as  many 
other  sinusoidal  wave  forms  (at  other  frequencies)  as  desired  to 
represent  an  offshore  spectrum,  and  the  results  added  together 
to  produce  a  time  history  of  the  observed  waves  at  a  fixed  point 
relative  to  the  moving  craft.  With  the  reference  offshore  wave¬ 


length  given  by  Aq 


where  u>m  is  the  frequency  of  maximum 


spectral  energy,  the  wave  pattern  changes  in  the  region  near  that 
defined  in  Equation  (25)  from  an  oscillatory  type  of  wave  to  a 
solitary  wave  form.  In  order  to  insure  continuity,  the  solitary 
wave  crest  follows  the  trough  of  the  last  oscillatory  wave  form 
in  this  "transition"  region,  with  a  continuous  numerical  transi¬ 
tion  weighting  of  the  oscillatory-type  waves  and  the  complete 
solitary  wave  form.  The  height  of  the  solitary  wave  form  in  this 
initial  transition  region  is  assumed  to  be  equal  to  the  signifi¬ 
cant  wave  height  of  the  original  offshore  wa'.es.  A  discussion 
of  the  computational  technique  for  achieving  this  model  of  a 
solitary  wave  i3  presented  in  a  later  section  of  the  report. 


For  the  present  problem,  the  wave  field  is  evaluated  in 
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terms  of  the  inertial  reference  coordinate  xQ,  which  is  measured 

from  an  initial  origin  at  distance  from  shore,  where  the  x  - 

axis  is  perpendicular  to  the  beach.  The  solitary  wave  propagates 
according  to 


n 


\ 

\ 

X  + 

4 

u  dt  -  ct  V 

(39) 


relative  to  the  ship,  for  a  following  sea  case,  with  x  measured 
relative  to  the  craft  CG.  The  height  change  follows  the  relation 
given  by  Equation  (27) ,  with  the  depth  variations  as  functions  of 
time  obtained  from  Equation  (37).  This  continues  until  the  limit 
condition  for  breaking  given  by  Equation  (26)  is  reached,  after 
which  the  only  effect  remaining  (after  breaking)  is  the  water 
elevation  equal  to  the  run-up  R,  which  can  be  represented  as  an 
increased  surface  elevation,  with  a  longitudinal  velocity  given 
by  Equation  (33). 


For  the  head  sea  case,  associated  with  retraction  from 
the  beach,  the  basic  wave  geometry  variations  in  accordance  with 
distance  along  the  beach  from  the  offshore  reference  point  remain 
the  same.  The  major  change  that  occurs  is  in  the  argument  of 
the  sech  function  of  the  solitary  wave  representation,  which  becomes 


(40) 


The  wave  heights,  breaker  forms,  etc.  associated  with  the  wave 
system  for  this  head  sea  case  are  obtained  in  the  same  manner 
as  for  the  following  sea  case. 

The  establishment  of  the  run-up  elevation  is  made  at 
the  point  of  breaking,  and  this  must  be  synchronized  with  the 
solitary  waves  at  breaking  to  insure  continuity  and  physical 
realism.  This  effect  must  also  be  re-initiated  at  each  solitary 
wave  breaking  to  provide  the  required  periodicity.  Methods 
of  carrying  out  these  operations  on  a  digital  computer  are 
provided  in  a  later  section. 

5.  EQUATIONS  Of  MOTION 

The  equations  of  motion  are  written  with  respect  to 
a  moving  axis  system  whose  origin  is  at  the  ship  CG.  These 
equations  represent  the  motion  of  the  ship  with  six  degrees  of 
freedom,  and  it  is  assumed  that  there  is  a  vertical  center- 
line  of  symmetry,  the  ship  mass  and  CG  location  are  constant 
for  any  mode  of  operation,  and  that  products  of  inertia  are 
zero.  For  the  present  case  of  the  assault  boat,  the  equations 
are  given  by  the  following: 


17 


NAVTRADEVCEN  68-C-0152-1 


m(u-vr+wq)  =  X 

(41) 

m(vtru-pw)  =  Y 

(42) 

• 

m(w-qu+pv)  =  Z 

(43) 

I  p  +  qr (I  -I  )  =  K 
xxr  ^  z  y 

(44) 

Iyy4  +  PrUx-V  =  M 

(45) 

1zzi  +  pqUy-V  -  n 

(46) 

The  external  forces  and  moments  acting  on  the  boat  (indicated  by 
the  right  hand  sides  of  the  equations  above)  arise  from  hull  hydro- 
dynamic  and  hydrostatic  reactions;  propeller  effects;  rudder 
deflection;  control  action;  and  aerodynamic  reactions  due  to  wind. 
These  are  represented  as 

F  =  F(hull)  +  F (props)  +  F (rudders)  +  F (control) 

+  F  (waves)  +  F(wind)  (47) 

and 

M  =  M(hull)  +  M (props)  +  M( rudders)  +  M( control) 

+  M  (waves)  +  M(wind)  (48) 


in  general  form  for  the  forces  and  moments. 


The  ship  motions  are  referred  to  an  inertial  reference 
frame  in  order  to  determine  the  absolute  position  relative  to  the 
shoreline  and/or  an  operating  area  offshore.  The  origin  of  the 
inertial  reference  frame  is  chosen  to  be  located  at  an  offshore 
point  1,000  yards  from  the  beach,  with  the  xQ-axis  directed  normal 

to  the  shoreline,  which  is  then  parallel  to  the  yQ-axis.  The  main 

angle  indicating  orientation  of  the  ship  with  respect  to  the 
inertial  reference  frame  in  the  present  case  is  the  yaw  angle  \p, 
since  the  roll  and  pitch  angles  are  expected  to  be  small.  The 
inertial  axis  system  and  the  body  axes  are  related  by  (see  [16]) 
the  matrix  transformation 
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— 
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x 

B 

B 

B 

X 

o 

XX 

xy 

xz 

y 

B 

B 

B 

y 

Jo 

yx 

yy 

yz 

z 

B 

B 

B 

z 

o 

zx 

zy 

zz 

_  ___ 

■ 

— 

_ 

(49) 


where 
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B  =  cos  0  cos  <p  7  B  =  sin  4-  sin  0  cos  ip  -  cos  4>  sin  ip 
xx  xy 

B  =  cos  4>  sin  0  cos  ip  +  sin  4>  sin  ip  7  B  =  cos  0  sin  4> 
xz  yx 


B  =  sin  ip  sin  0  sin  ip  +  cos  4  cos  ip 


(50) 


E,  =  cos  4>  sin  0  sin  ip  -  sin  <j>  cos  41  ;  B  =  -  sin  0 
yz  zx 

B  =  sin  4>  cos  0  7  Bzz  =  cos  4  cos  0 

For  small  values  of  <j>  and  0,  with  large  yaw  angle  ip  possible,  and 
retaining  only  linear  terms ,  the  relation  becomes 


MM  — ” 

““  — 

—  — 

/ 

xo 

cos  ip  -sin  4/  9  cos  <p  +  <j>  sin  ip 

X 

sin  ip  cos  ip  0  sin  4»  -  4>  cos  4* 

y 

zo 

-0  <p  1 

z 

(51) 


In  the  present  problem  involving  a  ship  on  the  surface, 
the  heave  motion  is  also  small,  and  hence  analysis  of  coupled 
heave  and  pitch  can  be  refined  to  a  fixed  orientation  axis  system 
that  translates  with  the  body  but  does  not  rotate  with  the  body. 

In  that  case  the  z-axis  displacement  for  the  body  axis  system,  for 
the  fixed  orientation  system,  and  also  for  the  inertial  system 
(zQ-axis)  are  identical,  within  linear  theory,  for  the  present 

purposes.  The  only  relations  between  the  coordinate  frames  are 
those  concerned  with  horizontal  plane  displacements,  and  these 
navigation  equations  are  then 

x  =  x  cos  ip  -  y  sin  ip 

(52) 

yQ  =  y  cos  +  x  sin  ip 
when  the  craft  is  proceeding  toward  shore,  and 

x  =  y  sin  ip  -  x  cos  i|» 

°  (53) 

y  =  -  y  cos  ip  -  x  sin  4/ 
o 

when  the  craft  is  retracting  from  the  bench.  An  illustration  of 
the  relations  between  the  inertial  frame  ard  the  ship  axis,  when 
the  ship  is  in  either  of  these  respective  modes  of  operation, 
is  shown  in  Figure  6.  For  the  case  where  the  boat  is  retracting 
from  the  beach,  in  a  stern-first  attitude,  the  body  axis  system 
on  the  boat  is  changed.  In  that  case  the  x-axis  is  positive  from 
the  CG  to  the  boat  stern,  and  other  axis  orientations  are  similarly 
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changed.  The  yav,  angle  f ,  when  measured  from  the  inertial  frame, 
is  really  (180°+']»)  ,  with  changes  only  considered  relative  to  an 
imaginary  axis  parallel  (but  opposite)  to  the  xQ~axis,  as  indicated 

in  Figure  6. 

In  the  further  derivation  of  the  equations ,  it  is  assumed 
that  the  vertical  plane  motions  of  heave  and  pitch  have  negligible 
coupling  with  the  lateral  plane  motions  of  yaw,  sway,  and  roll. 
Similarly  the  effect  of  surge  (i.e.  forward  speed  variation)  is 
coupled  with  the  lateral  plane  motions,  as  well  as  with  the 
vertical  plane  heave  and  pitch  motions.  These  assumptions 
regarding  relations  between  the  variables  allow  simplification 
of  the  left  hand  sides  of  Equations  (41) -(46)  and  in  the  evalua¬ 
tion  of  the  various  hull  forces  and  moments.  The  various  forms  in 
tlie  equations  are  obtained  by  analysis  (see  [2])  of  the  model  test 
data  in  [1],  and  there  are  often  significant  differences  in  the 
equations  for  the  two  configurations  considered,  due  to  the 
reversal  of  the  bow  and  stern  aspect  of  the  boat  when  it  travels 
toward  the  beach  and  retracts  seaward  from  the  beach.  In  the 
resulting  equations  of  motion  the  configurations  will  be  identi¬ 
fied  as  heavy  and  light,  with  the  understanding  that  the  light 
configuration  applies  to  the  retraction  stern-first  mode  of 
operation. 

On  the  basis  of  the  preceding  discussion,  the  equations 
of  motion  are  simplified  to 


m(u-vr)  =  X 

(54) 

m(v+ru)  =  Y 

(55) 

mz"  =  Z 

(56) 

I  p  -  K 
xx^ 

(57) 

I  q  =  M 

yy^ 

(58) 

I  r  =  N 
zz 

(59) 

In  order  to  complete  the  general  representation  of  the  equations 
of  motion,  the  following  kinematic  relations  are  included t 


_  _  d$  _  i 

p  -  at  "  ♦ 

de  • 
q  =  at  =  e 

d  if;  • 

r  =  at  =  * 


(60) 

(61) 

(62) 


where  <p  is  the  roll  angle  and  0  the  pitch  angle. 
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6.  HULL  FORCES  AND  MOMENTS 

The  forces  acting  on  the  hull  arise  from  both  hydro- 
dynamic  and  hydrostatic  effects,  where  the  hydrostatic  effects 
are  only  included  in  the  heave,  pitch,  and  roll  degrees  of 
freedom,  and  manifest  themselves  as  displacement-type  terms. 

The  hydrodynamic  forces  arise  due  to  inertial,  circulatory, 
and  viscous  effects  due  to  interaction  between  the  hull  and 
the  fluid.  The  forces  are  expressed  as  functions  of  the 
accelerations  and  velocities  of  the  hull,  and  for  the  motions 
of  surge,  sway  and  yaw  are  generally  represented  in  a  Taylor 
series  expansion.  A  number  of  terms  are  expected  to  be  small 
on  the  basis  of  symmetry,  and  hence  they  are  neglected,  and  the 
representation  is  primarily  dominated  by  linear  terms  in  the 
equations.  This  is  due  to  the  fact  that  only  small  angular 
changes  will  be  made  by  the  coxswain  during  the  course  of  his 
motion  toward  and  away  from  the  beach,  since  he  will  always  aim 
to  keep  the  boat  normal  to  the  wave  crests,  and  since  extensive 
studies,  e.g.  [3] ,  have  shown  that  the  motions  of  heave  and 
pitch  in  waves  are  primarily  linear. 

The  hull  forces  are  represented  in  terms  of  a  varying 
forward  speed,  since  that  quantity  will  be  controlled  during 
the  course  of  different  maneuvers.  As  a  result  certain  nonlinear 
terms  involving  forward  speed  will  appear  in  the  equations,  and 
these  are  primarily  parametric  in  nature.  The  hydrodynamic 
forces  due  to  hull  motions  are  considered  to  be  independent 
of  propeller  and  rudder  actions,  and  any  effect  on  net  forces 
is  taken  into  account  when  considering  those  particular  elements. 
Thus  the  axial  hull  force  terms  are  limited  in  this  section, 
since  the  data  analysis  [2]  of  the  test  results  in  [1]  con¬ 
sidered  the  net  effect  of  thrust  and  resistance,  ’-dilch  is 
treated  in  another  section. 

The  equations  for  hull  forces  for  the  surge,  sway  and 
yaw  modes  are  generally  expressed  as  factions  of  the  motion 
parameters  representing  the  velocities  and  accelerations  in 
these  modes  of  motion.  A  Taylor  series  expansion  is  carried 
out  up  to  second  order  in  order  to  include  nonlinear  effects, 
and  a  typical  example  for  the  case  of  the  lateral  hull  force 
is  shown  by 

Y (u, v,r,u, v,r)  =  Ye  +  YuAu  +  Yvv  +  Yrr  +  Y^v  +  Y^u  +  Y^r 

+  YvuvAu  +  YrurAu  (63) 

where  the  subscripts  represent  partial  derivative  operations, 
and  the  initial  equilibrium  condition  (denoted  by  subscript  e) 
is  straight  ahead  motion  at  constant  speed.  Since  there  are 
no  effects  on  the  Y- force  due  to  axial  acceleration  u;  no 
lateral  force  exists  in  equilibrium;  and  the  ship  is  almost 
symmetric  fore-and-aft,  a  number  of  coefficients  are  zero  (or 
negligible)  ,  i.e.  Y  -  Y  =  Y^  =  Y^,  ~  0.  With  the  forward  speed 
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being  a  total  variable,  no  perturbation  effects  due  to  speed 
changes  alone  exist,  and  since  (see  [16]) 


*v  -  I  l2  u  Xv 

Yr  =  |  L>  u  YJ. 


(64) 


where  Y'  and  Y'.  are  measured  nondimensional  coefficients,  all 

effects  of  forward  speed  and  the  variables  v  and  r  are  incor¬ 
porated  therein,  without  the  necessity  of  inclusion  of  any 
second  order  terms  involving  Au  (forward  speed  perturbation) .  As 
a  result  the  lateral  hull  force  is  represented  by 


Yhull  " 


Y  *  v  +  Y.uv  +  Y„ur 
v  1  2 


where 


(65) 


Y 


=  %  L2  Y1 
12  v 


(66) 


In  a  similar  manner  the  hull  yaw  moment  is  represented  by 


N 


hull 


=  N*r  +  N,ur  +  N_r  r  +  N_uv 


(67) 


where  the  rjrj  term  in  the  yaw  moment  equation  is  included  to 
allow  for  yaw  damping  effects  at  low  forward  speed,  when  turning 
around  during  retraction,  or  for  hull  forces  to  resist  propeller 
thrust  control  actions  when  beached.  The  form  of  the  hull  reactions 
Yhun  and  in  Equations  (65)  and  (67)  is  valid  for  both  the 

heavy  configuration  and  the  light  configuration,  but  different 
values  of  the  coefficients  will  generally  be  present  for  each 
configuration. 

The  axial  hull  force  expression,  neglecting  any  resistance 
terms  associated  with  purely  axial  motion  (which  are  included 
together  with  net  propulsive  forces)  ,  is  given  by 

Xhuii  =  X-u  +  Xjulvi  +  X2vr  (68) 

for  the  heavy  configuration,  and 

*hull  *  Xi“  +  X2vr  +  V*  (69) 

for  the  light  configuration.  The  major  difference  between  the 
two  expressions  is  the  form  of  the  axial  force  variation  due  to 
sideslip  angle,  i.e.  the  u!v|  ana  v?  terms,  which  occurs  due  to  the 
manner  of  mathematical  fit  to  the  available  experimental  data 
(see  [2]).  The  X^  term  is  the  longitudinal  added  mass,  and  the 

X2vr  term  occurs  as  a  hydrodynamic  inertial  reaction,  as  indicated 

by  potential  flow  dynamic  equations  in  [17]. 
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For  roll  motion  the  hull  equations  are  governed  by  forces 
due  to  roll  inertia,  damping,  and  hydrostatic  restoring  effects. 
Coupling  is  known  to  exist  with  yaw  and  sway,  with  the  main  effect 
due  to  sway,  while  roll  effects  (due  to  the  hull)  on  yaw  and  sway 
are  generally  negligible.  A  simplified  model  of  hull  roll  moments, 
originally  derived  in  [4] ,  is  used  herein  and  allows  for  the  coupling 
of  sway  in  modifying  the  resulting  coefficients  in  the  roll  equa¬ 
tion.  Thus  the  hull  roll  moment  representation  is  given  by 

Khun  =  K^’+  V  ‘  w)GMi*  '70) 

where  W  is  the  total  weight  (displacement)  and  | GM j  is  the 
metacentric  height. 

When  analyzing  roll,  and  also  for  heave  and  pitch,  which 
are  the  primary  wave-influenced  motions  due  to  their  hydrostatic 
restoring  terms,  the  equations  are  usually  written  with  respect 
to  an  axis  system  that  translates  with  the  body,  but  with  a 
fixed  orientation,  i.e.  they  do  not  rotate  with  the  body.  Since 
these  motions  are  linear,  and  do  not  couple  with  the  other  motions 
in  a  nonlinear  manner,  they  can  be  treated  in  their  usual  manner 
without  any  loss  of  generality  in  the  complete  six  degree  of 
freedom  equations.  In  the  derivation  of  the  equations  for  heave 
and  pitch,  the  boat  is  assumed  to  be  a  symmetrical  rectangular 
block  of  beam  B,  length  21,  and  draft  d,  with  the  CG  (origin  of 
coordinates)  assumed  to  be  at  the  geometric  center  (axially) 
and  the  value  of  the  draft  selected  as  an  average  that  provides 
the  correct  total  displacement  (d  =  2.7  ft.  for  the  heavy 
condition  and  d  =  1.4  ft.  for  the  light  load).  This  procedure 
is  sufficiently  accurate  for  representing • the  hull  forces  and 
moments,  allowing  simplification  due  to  symmetry  for  these  terms 
as  well  as  for  the  wave-induced  forces  (to  be  discussed  in  a 
later  section) . 

With  the  positive  direction  of  z  downward,  the  deriva¬ 
tion  of  the  local  inertial  hydrodynamic  vertical  force  at  a  section 
is  obtained  from  the  simplified  result  of  slender-body  theory 
[18] ,  which  states  that  the  local  force  on  any  section  is  equal 
to  the  negative  time  rate  of  change  of  fluid  momentum.  This  is 
expressed  by 


where  the  time  derivative  operation  is  defined  by 


(71) 


D_  _  3 _ 

Dt  at 


(72) 


A^3  is  the  vertical  added  mass  (two-dimensional)  of  the  cross- 
section,  and  w^  is  the  section  vertical  velocity  (along  z- 
direction,  i.e.  positive  down)  given  by 
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w 


b 


D 

Dt 


(z-50) 


Z-C6+U6 


(73) 


In  the  above  equations  the  coordinate  ?  is  a  "dummy”  variable 
along  the  longitudinal  body  coordinate  x,  and  coincident  with  it. 
Since  the  boat  is  assumed  to  be  a  uniform  rectangular  block, 

dA' 

the  value  of  A^  is  constant  at  every  section  and  g—  -  =  0  except 

at  the  ends,  where  it  has  6-function  form.  Expanding  the  operation 
of  Equation  (71)  and  integrating  over  the  hull  yields 


Z. 

l 


*£ 

A33 

. 

-£ 


u 


>£ 

A^  ® 

. 

— £ 


(74) 


£ 


where  the  symmetry  condition 


a!^  £;  dc  =  0  is  invoked. 


-£ 


The  damping  force  in  heave  arises  due  to  energy  dissipa¬ 
tion  associated  with  wave  generation  by  the  ship  motions  on  the 
free  surface.  The  ratio  of  the  amplitude  of  the  heave-generated 
two-dimensional  waves  to  the  heave  motion  amplitude  of  a  section 

is  denoted  as  S'  ,  and  by  equating  the  energy  balance  between  ship 
z 

work  and  outgoing  wave  energy,  the  vertical  damping  force  per  unit 
vertical  velocity  of  the  section  is  expressed  as 


pg2\2 


(75) 


where  w  is  the  circular  frequency  of  oscillation  (usually 

frequency  of  encounter  for  ship  motions) .  The  vertical  damping 
force  at  each  section  is  then 


-  N;2(z-?e+ue)  ,  (76) 

and  the  total  vertical  damping  force  is  determined  by  integrating 
over  the  ship  length,  leading  to 


Z 


d 


n;z  d< 


z  -  u 


N'  dc 
zz 


-x 


-£ 


(77) 


and  applying  symmetry,  since  N'  is  constant  for  each  section. 

z  z 
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The  hydrostatic  restoring  force  due  to  static  vertical 
displacement  (on  a  calm  free  surface}  of  a  ship  section  is 

dZ 

a~  =  -  PgB*(z-5e)  ,  (78) 


leading  to  a  total  force  given  by 


Z.  =  -  pg 


B*  d£  *  z  =  -  2pg£B  •  z 


-Z 


since  B*  =  3  is  constant  for  each  section. 


(79) 


The  pitch  moment  terms  corresponding  to  all  the  force 
terms  derived  in  the  foregoing  are  easily  obtained  by  use  of  the 
relation 


dM  _  dZ 

d£  "  "  ^  d£ 


(80) 


followed  by  the  required  integrations  over  the  hull.  This  leads 
to 


M. 


l 


(Z 

fl33  52  •  s  +  u 

A’ 3  d£  »  z  +  u2 

« 

4 

~Z 


-Z 


A33  ^ 


-Z 


(81) 
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(82) 


(83) 


The  added  mass  and  damping  coefficients,  A'  ^  and  N'  , 

33  zz 

appearing  in  the  preceding  equations  are  functions  of  the 
geometry  of  the  ship  section  (which  is  assumed  constant  in 
the  present  case  for  treating  heave  and  pitch)  and  are  also 
functions  of  frequency.  They  are  determined  from  solution  of  a 
hydrodynamic  free  surface  boundary  value  problem  for  a  class  of 
representative  ship  sections,  as  shown  by  the  results  in  [19]. 

In  order  to  allow  a  simple  representation  of  the  equations,  an 
approximation  is  made  in  selecting  constant  values  for  these 
terms  which  are  independent  of  frequency,  with  these  constants 
selected  in  a  manner  so  that  the  basic  dynamic  characteristics  of 
this  system  are  properly  exhibited.  Since  the  hydrostatic 
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restoring  terms  are  constant,  the  selection  of  the  added  mass  is 
such  that  the  natural  periods  in  heave  and  pitch  are  close  to 
those  observed  experimentally.  Similarly  the  damping  values  are 
chosen  at  values  appropriate  to  the  range  of  frequencies  that 
are  anticipated  for  the  assault  boat  when  it  approaches  the 
beach  or  when  retracting.  As  a  result  of  the  different  dis¬ 
placements  associated  with  these  two  operations,  and  the 
expected  range  of.  speeds  in  each  case,  a  different  set  of 
constant  coefficients  is  determined  for  the  two  ship  configura¬ 
tions.  An  outline  of  the  basis  for  selection  of  these  values  is 
provided  in  [2] .  Similar  reasoning  is  applied  to  the  determina¬ 
tion  of  the  a 'ied  roll  inertia,  while  the  roll  damping  is 
determined  from  combined  theory  and  empirical  data.  Reference 
[2]  provides  the  detailed  information  on  the  determination  of 
those  parameters  in  accordance  with  this  procedure. 

In  the  course  of  the  hydrodynamic  force  representations 
developed  above,  certain  assumptions  have  been  made  that  were 
not  delineated  in  detail,  and  they  are  described  by  the 
following.  Hydrodynamic  forces  due  to  inertia  asymmetries,  which 
were  explicitly  neglected  in  the  case  of  heave  and  pitch,  have 
also  been  neglected  in  the  representation  of  the  lateral  force 
ajid  yaw  moment.  In  particular  the  stability  derivative  terms 
Y*  and  N*  (in  conformity  with  SNAME  notation  [16])  have  been 

assumed  to  be  zero  and  have  been  deleted  from  the  listings  in 
Equations  (65)  and  (67) .  Another  assumption  made  is  the  lack 
of  any  influence  of  the  varying  depth,  along  the  sloping  beach, 
on  the  values  of  the  coefficients  included  in  the  previous 
equations.  Measurements  made  in  [20]  have  indicated  only  small 
variations  in  lateral  motion  hydrodynamic  coefficients  (yaw  and 
sway) ,  within  the  range  of  variation  considered  significant 
in  comparison  with  other  contributory  factors  such  as  frequency 
and  forward  speed.  Thus  it  is  assumed  herein  that  constant 
values  exist  for  these  coefficients,  and  that  the  varying  shallow 
water  influence  is  negligible.  Since  average  constant  coeffi¬ 
cients  have  been  selected  for  the  force  terms  for  heave,  pitch, 
and  roll,  a  similar  neglect  of  such  effects  is  inherent  in  the 
equations.  The  force  expression  that  is  known  to  vary  signi¬ 
ficantly  as  the  depth  changes  is  the  axial  force,  and  data  perti¬ 
nent  to  that  variation  was  obtained  in  [1]  and  analyzed  for  use 
in  the  present  report  in  [2],  This  data  will  be  discussed  while 
considering  propeller  effects,  since  the  net  axial  force  is  the 
quantity  of  important  for  the  equations  in  this  program. 

7.  PROPELLER  FORCES 

The  LCM  (6)  is  a  twin-screw  propeller  craft,  with  both 
propellers  being  right-handed  [5).  These  propellers  have  three 
blades,  with  2  ft.  diameters,  and  are  described  in  [6] .  The 
propeller  effects  acting  on  the  assault  boat  only  produce  an 
axial  force  (along  x~axis) ,  a  lateral  force  (along  y-axis) ,  and 
a  yawing  moment  (about  z-axis) .  The  present  section  treats 
the  case  where  both  propellers  are  operating  in  unison  to  propel 
the  craft,  and  any  considerations  of  differential  propeller 
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action  are  treated  in  the  section  concerned  with  control  effects, 
where  such  control  effects  are  due  to  the  action  of  the  coxswain 
via  his  control  of  the  propellers  (rudder  control  is  considered 
separately,  in  another  section).  From  the  tests  in  [1],  as 
analyzed  by  the  procedures  in  [2] ,  information  was  obtained  on  the 
variation  of  the  net  axial  force,  i.e.  T-R,  the  difference 
between  thrust  and  resistance,  for  a  range  of  ship  speeds  and 
propeller  rpm  values.  This  data  allowed  determination  of  the 
axial  force  as  a  function  of  ship  longitudinal  speed  and  propeller 
rotational  speed,  and  also  provided  data  for  determination  of  the 
self-propulsion  conditions,  i.e.  the  speed  corresponding  to  a 
particular  propeller  rotational  speed  for  twin  screw  operation. 

The  tests  carried  out  under  self -propulsion  conditions  then 
provided  information  on  the  side  force  and  yawing  moment  associ¬ 
ated  with  propeller  operation  under  these  conditions ,  from  which 
a  mathematical  form  was  fit  in  order  to  provide  expressions  for 
use  in  the  present  mathematical  model.  Self-propulsion  informa¬ 
tion  is  presented  for  the  two  configurations  in  Figures  7  and  8, 
where  the  rpm  data  is  for  propeller  rpm.  Since  there  is  a 
reduction  gear,  with  ratio  1.51,  employed  on  this  vehicle  [21], 
the  engine  rpm  is  increased  by  that  magnitude.  This  data  is 
applicable  only  to  deep  water  conditions  for  each  case. 

Since  the  propellers  are  right-handed,  a  starboard 
side  force  and  a  yawing  moment  to  port  are  obtained  for  forward 
operatic*,  (heavy  configuration)  and  these  terms  are  represented 
by 


Y (props )  =  ~7-13  u  +  5’74  u2  ~  °-75  u3  +  °-03  u4  (84) 

N (props)  =  103*13  u  "  114*2  u2  +  15-61  u3  ~  0.75  u4  (85) 

which  are  only  functions  of  the  forward  speed  (any  influence  of 
rpm  is  inherently  contained  within  the  speed  variation,  since 
there  is  a  unique  rpm  for  a  given  speed  under  normal  twin-screw 
operation) .  Similarly  the  side  force  and  yawing  moment  for  astern 
operation  (light  configuration)  are  given  by 


Y (props) 

=  6.91  u  +  0.58  u2 

(86) 

N (props) 

=  156.91  u  +  15.78  u2 

(87) 

where  the  axis  orientation  for  this  condition  is  indicated  in 
Figure  6.  It  is  this  effect  of  a  lateral  force  and  yaw  moment 
due  to  similar  rotational  directions  of  the  twin  propellers 
that  causes  coursekeeping  difficulty  for  this  craft. 

Ihe  net  longitudinal  force  acting  on  the  craft  as  a 
function  of  forward  speed  and  propeller  rpm  is  found  by- 
obtaining  a  "collapsed"  form  of  the  data,  in  the  same  manner 

T 

as  the  ordinary  thrust  coefficient  parameter  K_  =  — - — 

1  ,n2n'i 

K  1 
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for  a  propeller  is  represented.  The  resultant  form  is  somewhat 
complicated  due  to  the  wide  range  of  speeds  in  this  program, 
especially  for  the  forward  motion  of  the  heavy  configuration. 

The  force  expression  for  the  heavy  configuration  is  given  by 

X (props)  =  4'93  no2  +  °'82  no  "  °*06  u3  '  (88) 


where  nQ  =  operating  propeller  rotational  rate  (rps) .  This 

expression  is  found  by  means  of  fitting  the  data  with  a  mathe¬ 
matical  expression,  after  obtaining  the  collapsed  form  for  the 
range  of  operating  forward  speed  and  propeller  rotational  speed 
(see  [2]).  The  operating  propeller  rotational  rate  n  is 
defined  by  ° 

no  "  I  (Vns}  '  (89) 


which  is  the  mean  of  the  port  and  starboard  propeller  rates. 
When  both  propellers  are  operating  in  unison  at  the  same  rate 
the  value  of  nQ  is  proper,  as  required  in  the  present  section, 

and  any  effect  of  different  propeller  rates  (port  and  starboard) 
is  considered  when  treating  control  forces.  The  relation  in 
Equation  (88)  is  valid  for  the  case  of  motion  in  deep  water, 
where  no  influence  of  the  shallow  depth  is  present.  The  effect 
of  motion  along  a  beach  where  the  depth  is  decreasing  results 
in  an  increased  resistance,  which  is  found  (by  mathematical 
fit,  see  [2])  to  be  represented  in  the  form 


x<2> 

(props) 


0.08  un2  -  77.37 
o 


u2n  a 


(90) 


where  h  is  the  local  depth  at  the  ship  CG  and  a  is  the  beach 
slope.  The  sum  of  the  two  expressions  in  Equations  1 8 8 )  and  (90) 
is  the  axial  force  appropriate  to  the  operation  of  the  heavy 
configuration  when  approaching  a  shoreline. 

Similar  results  are  found  for  the  case  of  the  light 
configuration,  where  the  deep  water  axial  force  is  represented  by 

(1)  , 

X.  .  =  16.36  n2-  17.97  un  -  3.38  u2  (91) 

(props)  o  o 

and  the  effect  of  shallow  depth  is  given  by 


(■•) 

X  , 

(props; 


0.1:  jr. 

i  > 


u- n  . 

62.86  —  — 
i : 


(92) 
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8 .  RUDDER  FORCES 

The  effect  of  rudder  deflection  produces  a  side  force, 
yawing  moment  and  axial  force.  The  LCM  (6)  has  twin  rudders 
that  deflect  together  in  unison.  The  rudders  lie  in  the  wake  of 
the  propellers  when  the  assault  boat  moves  forward  in  the  heavy 
configuration,  while  they  are  ahead  of  the  propellers  in  the 
light  configuration  associated  with  the  retraction  operation. 
Thus  the  forces  and  moments  from  the  rudders  depend  upon  the 
rudder  angle,  the  longitudinal  ship  speed,  and  the  propeller 
rotational  speed,  with  less  dependence  on  the  propeller  speed 
for  the  light  configuration.  The  experimental  data  obtained 
in  [1]  was  analyzed  to  obtain  data  on  the  rudder  forces,  which 
were  found  to  depend  upon  forward  speed  wheri  represented  in  the 
form 


(rudders) 

-  n  i  L2»k6 

(93) 

(rudders) 

=  N1  |  L3u26 

(94) 

and  this  dependence  can  then  be  ascribed  to  the  effect  of  the 
propellers.  The  hydrodynamic  derivatives  Y'  and  N'  depend  upon 
forward  speed,  and  are  given  by  ' 

Y'6  =  0.00093  u  (95) 

and 

N’6  =  -0.00036  u  (96) 

for  the  heavy  condition,  and  by 

Y’6  =  0.00016  u  (97) 

and 

NJj  =  0.000078  u  (98) 

for  the  light  condition.  The  notations,  positive  direction  of 
rotation,  etc.  used  here  are  in  conformity  with  (16).  The  final 
representation  of  the  lateral  force  and  yaw  moment  for  rudder 
deflection  is  then  given  by 


Y 

(rudders) 

=  2,71  u35 

(99) 

and 

N  (rudders) 

=  -56.69  u 3 ^ 

(100) 

for  the 

heavy  configuration, 

and  by 

Y 

(rudders ) 

=  0.47  u3S 

(101) 

and 

^  (rudders ) 

=  12.28  u3 

'102) 
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for  the  light  configuration. 

Operational  experience  has  shown  that  the  rudders  are  not 
very  effective  in  controlling  this  craft  when  it  is  moving 
astern  (retraction  conditions) ,  and  hence  the  rudder  forces  are 
not  significant  in  comparison  to  the  other  forces  acting  upon 
the  ship,  at  the  low  speed  associated  with  this  condition  [5], 

As  a  result,  extensive  control  action  is  obtained  by  means  of 
controlling  the  propeller  thrust  in  a  differential  manner,  and 
the  rudder  effects  are  only  additional  to  this  main  method  of 
control  for  this  operational  condition.  The  relations  provided 
in  this  section  are  assumed  to  be  valid  for  rudder  angles  in  the 
range  up  to  ±30°,  where  linear  behavior  is  still  a  valid  approxi¬ 
mation  for  the  lateral  force  and  moment  (see  [2]). 

Some  consideration  was  also  given  to  an  alternate  repre¬ 
sentation  of  the  rudder  forces  in  a  manner  more  dependent  upon 
propeller  rotational  rate,  rather  than  forward  speed  alone.  While 
such  a  relation  was  obtained  (see  [2]  for  details) ,  the  present 
form  in  Equations  (99)- (102)  is  retained  herein  since  it  is 
anticipated  that  the  rudders  will  primarily  be  used  during  ahead 
motion  of  the  heavy  configuration  at  appreciable  forward  speeds. 
The  representation  in  terms  of  propeller  rotational  speed  (for 
the  heavy  configuration)  is  useful  in  low  speed  operations  where 
propeller-rudder  interaction  could  produce  turning  action  on  the 
craft,  and  that  form  is  presented  in  [2]  in  case  a  particular 
situation  arises  where  such  an  expression  would  be  more  useful 
in  certain  simulation  exercises. 


The  rudder  deflection  will  also  produce  an  additional 
resistance,  resulting  in  an  axial  force  proportional  to  <$2. 
This  force  is  expected  to  be  small,  and  hence  an  approximate 
representation  is  given  by 


X (rudders)  "  "2,S  uZfi2 

which  is  obtained  by  analytical  estimation  in  [2] . 


(103) 


Due  to  the  small  size  of  the  rudders,  and  the  direct 
manual  control  exercised  by  the  coxswain,  very  little  time  lag 
exists  in  actuating  the  rudder  to  a  commanded  value.  A  simple 
analysis,  based  on  considerations  ordinarily  used  in  treating 
such  effects  (22]  ,  shows  a  time  lag  of  less  than  1  second 
(-  .3  sec.  for  a  10  knot  forward  speed)  so  that  such  effects 
are  neglected  in  this  report. 


9 .  CONTROL  FORCES 


As  indicated  previously,  control  forces  are  obtained 
by  use  of  differential  thrust  of  the  propellers,  under  command 
of  the  coxswain.  While  it  is  known  that  the  lateral  force  and 
yaw  moment  effects  of  a  single  propeller  are  not  the  same,  when 
considering  the  separate  forces  exerted  by  either  the  port  or 
starboard  propeller  acting  alon  (see  the  results  in  [1]  and 
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discussion  in  [2]),  the  separate  thrust  force  effects  of  each 
propeller  can  be  considered  to  act  independently  of  each  other, 
as  a  first  assumption.  Thus  it  is  possible  to  represent  these 
effects  using  data  on  propeller  thrust  characteristics  alone  as 
obtained  in  [1] ,  with  the  assumption  that  equal  thrusts  are 
developed  by  each  propeller,  at  the  same  forward  speed  and 
rotational  speed,  when  operating  together  in  a  twin-screw 
installation. 


Measurements  of  the  propeller  thrust,  under  twin-screw 
action,  were  obtained  in  [1]  and  analyzed  in  [2] .  It  was  found 
that  the  total  thrust  under  these  conditions,  in  ahead  operation 
for  the  heavy  configuration,  can  be  represented  by 

T  =  13.44  n  2  ,  (104) 

o  ' 

Similarly  the  total  thrust  under  twin-screw  operation,  for  the 
light  configuration  moving  stern  first,  is  given  by 

T  =  0.95  nQ3  .  (105) 

These  results  are  appropriate  to  the  "normal"  operating  condition 
of  the  boat  when  moving  ahead  with  twin-screw  operation  in  unison, 
and  it  is  assumed  that  half  of  the  total  thrust  is  developed  by 
each  propeller. 


When  considering  differential  thrust  action,  a  rotational 
speed  different  from  the  normal  operating  value  referred  to  above 
is  selected  by  the  coxswain,  and  the  change  in  thrust  for  the  port 
or  starboard  propeller  (based  upon  the  reference  to  these  locations 
when  considering  the  boat  in  its  normal  bow-first  ahead  operation) 
can  be  represented  by 


where 


4Tp,s  -  I  Ks>  -  T(no>]  (106) 

n  =  rps  of  port  or  starboard  propeller, 
as  selected  by  the  coxswain  for 
differential  action  (107) 


and  n  is  defined  in  Equation  (89).  The  effect  of  this  altered 
thrust  generated  by  the  propellers  is  a  yawing  moment  due  to  the 
location  of  each  propeller  at  a  distance  of  2.5  ft.  from  the 
longitudinal  centerline  given  by 


N (control ’  1-25  [TV  -  T(no>] 
for  the  action  of  the  port  propeller,  and 

N  (control)  s  ■  -  i’25  [T("s>  '  T<"o> 


(108) 


(109) 
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for  the  action  of  the  starboard  propeller,  when  applied  to  the 
heavy  configuration  in  normal  forward  motion.  The  resulting  total 
control  yaw  moment  is  the  sum  of  these  expressions,  given  by 

"(control)  =  i'25  [fV  -  T‘V]  (110> 

which  only  depends  upon  the  values  of  the  port  and  starboard 
propeller  rotational  speeds.  Similarly,  for  the  light  configura¬ 
tion  appropriate  to  retraction  conditions,  the  yawing  moment  due 
to  differential  propeller  control  is  given  by 

“(control)  -  !-25  E<"s>  ‘  T<"p>]  <1U> 

The  effect  of  any  change  in  either  propeller  rotational  speed  on 
the  axial  force  is  manifested  by  the  use  of  the  resulting  value 
of  the  operating  propeller  speed  nQ  (defined  in  Equation  (89)), 

in  the  expressions  given  by  Equations  (88),  (90) -(92). 

Since  the  relatively  small  engines  of  the  LCM  (6)  are 
under  the  direct  control  of  the  coxswain  by  means  of  hand  throttles, 
it  is  expected  that  rapid  response  to  commanded  settings  of  the 
engine  will  occur.  Thus  the  resulting  propeller  rotational  speed 
will  respond  rapidly,  and  for  the  purposes  of  this  report  the  lags 
will  be  considered  to  be  negligible.  Since  the  propeller  force 
effects  due  to  changing  propeller  speed  will  produce  relatively 
small  accelerations  on  the  boat,  any  effects  due  to  neglecting 
the  small  time  lag  are  still  expected  to  result  in  smooth  motion 
outputs. 


10.  WAVE  EXCITING  FORCES  AND  MOMENTS 

The  wave-induced  forces  and  moments  acting  on  the  assault 
boat  are  made  up  of  hydrostatic  terms  due  to  periodic  buoyancy 
alterations  as  the  waves  progress  past  the  ship  hull,  and  hydro- 
dynamic  forces  due  to  inertial  effects  of  the  wave  velocities  and 
accelerations  acting  at  each  ship  section.  The  buoyancy  effect 
for  the  vertical  force  is  represented  by 

(i) 

(waves)  „  _pgB,„u,t)  (112) 


at  each  section,  and  the  inertial  contribution  to  the  wave 
exciting  forces  may  generally  be  represented  by 


where  A’.. 


n 

sectional 


dF 


D 


_ (waves)  =  D,  ,  }  4  dS  __ 

df.  Dt  lAii  o:'  4  pb  Dt 


is  a  particular  local  added  mass  term,  S 
area,  and  Vq  is  a  component  of  the  wave 


V  (113) 

o . 

i 

=  B*d  is  the 
orbital 
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velocity  evaluated  at  the  undisturbed  surface  along  the  ship 
vertical  centerline.  The  waves  are  generally  assumed  to  be 
propagating  in  a  direction  relative  to  the  ship  hull  defined  by 
the  angle  y,  where  y  is  the  angle  between  the  ship  x-axis  and  the 
normal  to  the  wave  crests.  For  the  present  case  the  angle  y  is 
defined  by 

Y  =  i|>  -  8  (114) 


where  8  is  the  angle  made  by  the  normal  to  the  wave  crests,  measured 
relative  to  the  x^-axis,  which  changes  as  the  shoreline  is 

approached  due  to  refraction  effects. 


The  surface  wave  elevation  relative  to  the  ship  is 
represented  by 


■t 

n  •-  a  sin 

k  £  cos  y  +  k  y  sin  y  -  wt  k  cos  y 

u  dt 

* 

o 

(115) 


in  terms  of  the  dummy  variable  £  along  the  ship  x-axis,  where  the 
amplitude  a  changes  as  the  boat  approaches  the  shore  in  accordance 


t 


with  Equation  (5)  (a  term  involving 


v  dt  has  been  neglected  in 


1  o 

Equation  (115)  as  its  influence  is  not  significant  for  wave  forces) . 
The  wave  orbital  velocities  along  the  x,  y,  and  z-axes  are  defined 
as 


uq  =  aw  coth  k.h  cos  y  sinl^  k£  cos  y  -  wt  +  k  cos  y 


u  dt  (116) 


vQ  =  aw  coth  kh  sin  y  sin 


o 

rt 


k£  cos  y  -  wt  +  k  cos  y 


udt  ] 


(117) 


w  =  aw  cos 
o 


! 

ft 

k£  cos  y  -  wt  +  k  COS  y 

u  dt 

i 

o  < 

(118) 


where  they  are  evaluated  at  the  mean  free  surface  level  along  the 
ship  centerline,  y  =  0.  The  orbital  accelerations,  which  are 
required  in  evaluating  the  wave  forces  according  to  Equation  (113) 
are  obtained  by  carrying  out  the  operation 


_  3 _ 

Dt  at 


u 


*  r 


(119) 


and  they  are  given  by 
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« 

uo  =  -agk  cos  y  cos  (k  £  cos  y  - 

f  (t) ) 

(120) 

* 

v  =  -agk  sin  y  cos  (k  £  cos  y  - 
o 

f  (t) ) 

(121) 

wq  =  aw2  sin  (k £  cos  y  -  f (t) ) 

(122) 

rt 

where  f (t)  =  wt  -  k  cos  y  u  dt  . 

•  o 

(123) 

The  wave  force  in  the  z-direction  due  to  hydrostatic 
effects,  as  indicated  by  Equation  (112),  is  given  by 


(i) 

L  (waves) 


-pgaB 


(l 

sin  (k  £  cos  y  -  f  (t) )  d£ 


-2pgaB 


cos  (k  £  cos  y)  d£  *  sin  f (t) 


*  o 

=  -2pgae  5i£LJ|A£2=ja  sin  £(t)  ,  (124) 

where  symmetry  considerations  have  been  applied.  The  inertial 
hydrodynamic  contribution  to  this  force,  according  to  Equation 
(113) ,  is  then  found  to  be 


Z(2) 

(waves) 


2  (pBd  +  A'j3) 


?  sin  (k  a  cos  y ) 

aw  - 1 - - - 1 — 

k  cos  y 


+  2aiouA'33  sin  (k  £  cos  y) 


sin  f(t) 


(125) 


where  the  o-function  behavior  of  the  added  mass  at  the  ends  is 
accounted  for  in  the  last  term  in  this  expression.  With  the 
pitch  moment  related  to  the  vertical  force  in  accordance  with 
Equation  (80),  it  can  be  shown  that  the  results  for  the  total 
wave-induced  pitch  moment  are 
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M  (waves )  ~  < 


2pgaB 


sin  (k  £  cos  y)  _  £  cos  (k  £  cos  y) 
L  (k  cos  y) 2  k  cos  y 


-  2{pBd  +  A'^)  aw2 


sin  (k^.  cos  y)  _  £  cos  (k  £  cos  y) 
(k  cos  y) 2  k  cos  y 


-  2aiuuA22  £  cos  (k  si  cos  y)  >  cos  f(t) 


(126) 


The  axial  and  lateral  forces  due  to  waves,  in  accordance 
with  Equations  (113),  (116),  and  (117),  are  then  given  by 

X (waves)  =  “  ^pgaBd  sin  (k  si  cos  y)  cos  f  (t)  (127) 


and 


Y 

(waves ) 


2  (pBd  +  A'22)  ag 


sin  (k  £  cos  y) 
cos  y 


+  2aa)  coth  kh  u  A'22  sin 


(k  £  cos 


Y> 


•  sin  y  cos  f (t) 


(128) 


where  A22  is  the  lateral  sectional  added  mass  of  the  rectangular 

block  representing  the  craft.  From  Equations  (113)  and  (128)  the 
yawing  moment  due  to  waves,  defined  by 


dN ,  .  dY ,  . 

(waves)  _  (waves) 

a? 


(129) 


is  then  given  by 


N  .  .  =  -<2(pBd  +  A'  )  — 

(waves)  ]  22  cos  y 


sin  (k  £  cos  y)  ,,  „  , 

- 1  £  cos  (k  £  cos  y  ) 


+  2aw  coth  hk  u  A'22  £  cos  (k  £  cos  y) 


sin  y  •  sin  f (t) 


(130) 
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With  the  simplified  form  of  roll  moments  due  to  the  hull 
given  by  Equation  (70) ,  the  roll  moment  due  to  waves  can  be 
represented  by 


K 


(waves) 


=  pgak  sin  y  | GM 


S  cos  (k  5  cos  Y  -  f (t)  d£ 


-Jc 


=  2pga  tan  y  | GM |  Bd  sin  (k  i  cos  y)  cos  f  (t)  ,  (131) 


which  is  based  on  a  simplification  of  the  wave-induced  roll 
moment,  as  shown  in  [4] .  This  result  is  sufficiently  accurate 
for  use  in  the  present  study  where  roll  is  not  a  predominant 
motion  of  concern. 


In  all  of  the  preceding  results,  the  application  has 
been  made  to  oscillatory  waves  propagating  toward  a  shoreline, 
and  the  magnitude  of  the  associated  parameters  such  as  the  wave 
amplitude,  wavelength,  etc.  will  vary  in  accordance  with  the 
position  of  the  craft  along  the  sloping  beach  (i.e.  as  a  function 
of  the  local  depth  h) .  The  force  and  moment  magnitudes  depend 
upc'i  the  angle  y,  but  during  the  course  of  approaching  a  beach 
the  coxswain  often  steers  so  that  the  boat  tends  to  be  normal 
to  the  crests,  resulting  in  a  small  value  of  the  angle  y.  Some 
simplifications  may  be  made,  depending  upon  the  magnitudes  of  y 
expected  in  accordance  with  the  initial  heading  direction  of  waves 
far  offshore,  so  that  the  approximations  cos  y  =  1  and  sin  y  -  y 
can  be  applied  in  some  cases. 


The  results  above  are  in  a  form  that  is  primarily 
applicable  to  the  case  of  following  seas,  when  approaching  the 
beach.  Application  to  the  head  sea  case,  corresponding  to 
retraction  when  returning  seaward,  is  made  by  the  appropriate 
change  in  the  definition  of  the  function  f(t)  to 


f  (t) 


ft 


wt  +  k  cos 


u  dt 


(132) 


changing  its  sign  wherever  it  appears.  The  resulting  change  is  a 
change  in  sign  of  the  Z(waves)  and  N (waves)  exPressions'  and 

also  a  change  in  the  expression  for  the  surface  wave  elevation. 
These  changes  are  necessary  as  a  result  of  the  change  in  body 
axis  orientations  for  the  two  operating  conditions,  as  indicated 
in  Figure  6,  and  they  will  produce  consistent  results  for  the 
motions  with  proper  phase  relations  maintained. 


Another  important  wa\e  force  necessary  for  this  investi¬ 
gation  is  the  added  resistance  due  to  waves,  which  is  a  nonlinear 
phenomenon.  Considerations  have  only  been  given  to  determining 
the  mean  value  of  this  force,  and  no  simple  expression  exists 
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for  theoretical  evaluation  of  either  the  mean  value  or  the  total 
time  variation  of  this  term.  Since  prior  experimental  work  [23] 
has  indicated  a  proportionality  of  this  force  with  the  square  of 
the  wave  amplitude  for  regular  waves,  a  model  for  this  added 
resistance  is  formulated  as 

=  -  A  Jn  l  (133) 


where  the 


H 


(waves) 


symbol  represents  the  envelope  of  the  encountered 


wave  time  history  (at  the  ship  CG) .  The  envelope  of  a  function  is 
obtained  by  a  full  wave  rectifying  operation  (i.e.  absolute  value 
of  the  function) ,  followed  by  a  low-pass  filter,  and  the  resultant 
function  produced  by  the  operations  of  Equation  (133)  will  be 
slowly  varying,  with  some  higher  frequencies  (of  small  magnitude) 
superposed  on  the  main  modulating  curve.  Since  the  motion  in 
head  seas  is  known  to  produce  the  largest  added  resistance,  the 
coefficient  values  (i.e.  A  values)  in  Equation  (133)  should 
reflect  the  difference  between  head  and  following  sea  results 
for  the  two  configurations.  For  the  following  sea  case  (heavy 
configuration)  the  expression  for  this  force  is 

=  -700  |nj  (134) 


(waves) 


(134) 


and  for  the  light  configuration,  in  head  seas 

(2)  f  ]2 

X (waves )  *  -300  " 


it  is 


(135) 


A  discussion  of  the  choice  of  these  values  and  their  relation  to 
observed  craft  behavior  in  model  tests  is  given  in  [2] . 


Since  the  waves  change  form  as  they  progress  toward  the 
shore,  and  are  represented  by  solitary  waves  in  the  present 
approximate  procedure,  the  wave  forces  acting  on  the  ship  will 
also  change  their  form.  The  expressions  for  the  wave  forces 
due  to  a  solitary  wave  are  obtained  by  the  same  procedure  as  in 
[24] ,  together  with  other  hydrodynamic  effects  not  included  in 
[24].  The  wave  forces  derived  in  [24]  only  consider  forces  for 
determining  surge,  heave,  and  pitch  motions,  and  they  are  derived 
on  the  basis  of  the  Froude/Krilof f  hypothesis  that  neglects 
interaction  between  the  wave  and  the  ship  (i.e.  using  the 
pressures  existing  in  the  wave  field  without  allowing  any  modifi¬ 
cation  due  to  the  presence  of  the  ship) .  The  hydrodynamic  inertial 
effects  in  waves,  represented  by  Equation  (113) ,  provide  additional 
force  components  for  those  modes  and  also  allow  extension  of  the 
results  of  [24]  to  provide  wave  force  and  moment  expressions  for 
the  other  three  motions  not  considered  in  [24]. 


Using  the  results  presented  in  [11]  for  the  solitary 
waves,  based  upon  the  wave  form  expression  given  in  Equation  (39) 
for  following  seas,  the  water  particle  velocities  evaluated  at  the 
still  water  level  along  the  ship  centerline  are 
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(136) 


(138) 


alon.,  the  x,  y,  and  z-axes  of  the  ship,  where  c  is  defined  in 
Equation  (29).  These  results  are  generalizations  for  a  heading 
between  ship  and  wave,  for  a  moving  reference  point,  with,  the 
elevation  of  the  solitary  wave  (along  the  ship  centerline,  y  =  0) 
represented  by 


tanh 


V4h3  L 


X  COS  Y  -  Ct  +  COS  Y  U  dt 


n 


where 


(139) 


(140) 


The  hydrostatic  pressure  at  any  depth  in  a  solitary  wave 

region  is 

p  =  og(n  +zQ)  ,  (141) 

where  zq  is  the  depth  below  the  undisturbed  free  surface  level. 

The  time-dependent  part  of  the  pressure  is  pgn,  and  this  quantity 
is  used  tr  derive  the  wave  forces  in  [24].  Generalizing  the 
result  in  [24]  to  the  present  case  results  in 
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X (waves)  =  j  sech2  e  cos  y  -  g(t)J 

-  sech2  e  cos  y  + 

for  the  longitudinal  force  in  a  solitary  wave,  where 

-  -\  /3B~ 

^  4h3 

The  vertical  force  in  a  solitary  wave,  from  [24]  is  given  by 

/ 

2  (wives)  "  -  ^  '  tanh  e  cos  y  g(t)] 

+  tanh  e  |1  cos  y  +  g  (t  }^j\ 
and  the  pitch  moment  by 


(142) 


(143) 


(144) 


(waves) 


_  PggBfc.  |  tanh  e  |~i  cos  y  -  g(t)~j 
-  tanh  e  £.  cos  y  +  g  (t)  > 

"  r  r-  J 


+  log 
c2 


cosh  e 

i  cos  y  +  g(t)1 

L  J 

cosh  t 

a  cos  y  -  g(t) 

-  j  pgBd2  ^sech2  e  |”fc  cos  y  -  g  (t)j 


-  sech2  e 


cos  y  +  g  (tfj| 


(145) 


The  results  for  the  wave  forces  in  Equations  (144)  and 
(145)  are  due  to  hydrostatic  effects,  while  the  longitudinal 
force  is  also  interpreted  in  that  way.  Additional  terms  due 
to  hydrodynamic  inertial  effects  exist  for  the  vertical  force 
and  pitch  moment,  while  these  hydrodynamic  terms  are  the  only 
proper  means  for  determining  the  lateral  force  and  yawing 
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moment  due  to  waves.  The  roll  moment  is  found  to  be  primarily 
determined  by  hydrostatic  effects,  after  elimination  of  cross¬ 
coupling  hydrodynamic  terms  in  the  equations  of  motion  (see 
[4]),  and  the  influence  of  hydrodynamic  interaction  for  the 
longitudinal  wave  force  is  small  due  to  the  small  longitudinal 
added  mass  of  the  boat  relative  to  its  physical  mass. 

Assuming  that  the  relation  in  Equation  (113)  can  be 
applied  to  the  present  shallow  water  case  with  a  solitary  wave 
representation,  the  wave  force  components  due  to  this  hydro- 
dynamic  effect  are  given  below  as 


ZU)  =  /3gK 


H 
i  h 


(pBd  +  A‘  )c 

- rrr— - A'  u 

cos  y  33 


•  jsech2  e  fx,  cos  y  ~  9  (t/J 


l 


•  tanh  e[t  cos  y  ~  g(t)J  +  sech2  e  cos  y  +  g(t)J 

•  tanh  £  j~£  cos  y  +  g  (t) 


) 


1 


I 


(14S) 


M 


(2) 


-  -  /3gH-  |  "  c 


3/2  (pBd  +  A«33) 


COS  Y 


•  sech2  e  £  cos  y  “  g(t) 


•  tanh  e  £  cos  y  ~  g(t)J  “  Xsech2  e  |^£  cos  y  +  g(t)J 


[' 


•  tanh  e  £  cos  y  +  g(t) 
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COS  Y 
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sech^  e  £  cos  y  ~  g(t) 


■  2~cb¥  Y  sech2  e|^  cos  Y  +  g(tm  +  /IgK  ^J72  A’33£U 


sech2  c 


H 


(waves)  h 


cos  y  -  g  (t) 

|  •  tanh  e 

|7  cos  y  -  g(tfl 

cos  y  +  g  (t) 

♦  tanh  e 

j  £  COS  Y  +  g(t)J 

c  tan  y  (pBd 

+  A33}  ’  I 

sech2  e  |7  cos  y 

(147) 


,  H  , 
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( 

J  noo hi  <-  i 


■/gF  A^2  u  sin  y  • 

J  H 

*./  —  rr  i.Tl  —  C?  A  u2  ,  17  /%  r*  ..  J*  «*  /  1.  \  li 

—||{148) 


0  r*/~%G:  •./  —  rr  If  I  l  —  <?«  . 

~ - -  f  v  Xw,  .  c.  I  *  J  I  Kj 


40 


NAVTRADEVCEN  68-00152-1 


N 


H 


r 


(waves)  =  -  R-  °  tan  Y  (pBd  +  A^)  •  U  sech*  c 
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-  sech2  ;  cos  y  +  g(t)J  l 


(150) 


Since  the  solitary  wave  action  occurs  near  shore,  where  the  boat 
will  generally  be  controlled  to  keep  the  angle  y  small,  the 
approximations  cos  y  «  1,  sin  y  “  tan  y  -  y  can  be  made  to 
simplify  the  foregoing  expressions.  When  applying  these  results 
to  the  head  sea  retraction  operation,  the  solitary  wave  form  is 
changed  to 


V— 

— 

~l 

1 

X  COS  y  +  ct  +  COS  Y 

u  dt 

1 

V4h3 

__ 

'  > 

J 

(151) 


the  sign  of  u  is  changed  in  the  arguments  of  hyperbolic  functions, 
and  the  signs  of  Z(waves)  and  N(waves)  are  a^ao  changed,  iust  as 

in  the  case  of  the  treatment  of  oscillatory  waves  discussed 
previously. 


The  forces  and  moments  due  to  a  solitary  wave  decrease 
rapidly  for  large  values  of  time,  and  heno  ■  they  are  similar 
in  form  to  transient  impulse  effects.  Tnis  behavior  is  appropriate 
to  the  treatment  of  breaker  effects  physically,  and  the  short  time 
of  action  allows  sim^  e  solution  of  the  motion  responses.  The 
only  problem  associated  with  the  treatment  of  these  forces  is  the 
continuous  variation  that  must  be  computed  for  the  parameters 
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such  as  H,  h,  etc.  as  the  ship  moves,  and  the  relation  of  the 
wave  geometry  to  the  ship  location  at  every  instant.  Similarly 
it  is  necessary  to  represent  some  type  of  "transition"  of  the 
wave  forces  as  the  wave  field  changes  from  oscillatory  to 
solitary  wave  form,  and  this  requires  maintaining  a  continuity 
in  time  of  the  various  forces  acting  on  the  boat  during  the 
transition  period.  These  same  requirements  and  problems  exist 
for  the  case  of  following  seas  when  approaching  the  beach,  and 
also  head  seas  for  retraction,  with  the  time  and  space  reference 
relations  for  boat  and  wave  field  as  primary  considerations. 

A  discussion  of  computational  techniques  for  generating  the 
solitary  wave  field  and  the  forces  associated  with  it  is 
given  in  a  later  section  of  this  report. 

The  forces  due  to  the  run-up  water  elevation  following 
breaking  are  approximated  as  equivalent  to  those  associated  with 
a  sudden  elevation  of  the  craft  to  a  new  height,  together  with 
a  force  in  the  x-direction,  a  force  in  the  y-direction,  and  a 
yawing  moment.  The  longitudinal  velocity  field  in  the  run-up 
bore  produces  a  net  flow  relative  to  the  craft,  given  by  (u  +u) , 

which  is  interpreted  as  an  equivalent  current,  where  uR  is  defined 

in  Equation  (33) .  The  axial  force  produced  by  this  flow  is 
considered  to  be  equivalent  to  a  resistance  force,  which  is 
assumed  to  be  frictional,  and  which  can  act  to  propel  the  ship 
forward  in  a  following  sea,  as  long  as  the  ship  is  in  the  flow 
field,  and  it  adds  additional  resistance  in  a  head  sea  situation. 
The  method  of  synchronizing  the  run-up  flow  with  the  solitary 
waves,  and  for  evaluating  forces,  is  presented  in  a  later 
section  of  the  report. 

The  longitudinal  force  is  represented  by 

XR=|CfS  (Uj^-u)2  cos  Y  (152) 

for  the  following  sea  case,  where  is  the  friction  coefficient 
and  S  the  underwater  wetted  surface  area,  and  by 

XR  =  -  j  Cf  S  (uR+u) 2  cos  y  (153) 

in  the  head  sea  case.  When  the  craft  is  beached  and  has  no 
forward  speed,  the  same  expression  as  Equation  (147),  with  u  =  0, 
is  valid  and  represents  a  force  that  occurs  at  the  same  periodicity 
as  the  run-up  flow  field.  The  use  of  these  expressions  neglects 
any  short  term  impulse-like  action  due  to  the  initial  contact 
of  this  flow  field  with  the  craft  hull,  and  is  an  approximation 
in  that  sense. 

The  transverse  force  in  the  run-up  field  is  assumed  to 
be  the  same  as  that  due  to  a  current  acting  on  the  craft.  On  the 
basis  of  analysis  of  deep  water  model  tests  as  presented  in  [25] , 
this  lateral  force  is  given  by 
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Yr  =  .6  Ld  (uR±u)2  sin  Y  (154) 

where  the  +  sign  is  for  head  seas  and  the  -  sign  for  following 
seas.  The  appropriate  value  of  d  for  each  configuration,  and 
the  value  of  y  for  the  particular  case  of  either  head  or  following 
seas  (as  evaluated  during  the  simulation)  determines  the  value 
of  the  lateral  force  in  this  case.  The  results  are  also  valid 
for  the  situation  where  the  craft  is  beached,  with  u  =  0.  The 
yawing  moment,  also  found  from  [25]  ,  is  represented  by 

Nr  =  .035  L2  d  (uR±u)2  sin  2y  (155) 

with  the  same  interpretation  of  signs,  parameter  values,  etc. 
as  for  the  lateral  force.  In  all  of  the  preceding  expressions 
it  is  necessary  to  relate  the  run-up  flow  field  to  the  craft 
geometrically,  since  no  force  will  act  on  the  craft  until  both 
craft  and  run-up  flow  are  in  coincidence  and  the  flow  field  acts 
upon  the  craft.  Further  discussion  of  this  requirement  will  be 
presented  in  a  later  section. 

The  vertical  force  and  pitch  moment  due  to  the  run-up 
is  approximated  by  the  hydrostatic  effect  alone,  as  indicated  e.g. 
by  Equation  (112) ,  with  special  mathematical  functions  represent¬ 
ing  the  run-up  as  a  function  of  time  and  space.  With  the  run-up 
R  assumed  to  have  a  constant  magnitude  for  a  given  offshore  sea 
condition,  and  the  constant  ship  beam,  the  vertical  force  and 
pitch  moment  are  given  by 


Z 


R 


-pgBR 


I  fR(e,t)  d£ 
*  hull 


(156) 


and 


-  p  gBR 


e  fR(5/t)  d£ 


hull 


(157) 


where  f_(£,t)  is  a  special  function  describing  the  time  and  space 

X\ 

dependence  of  the  run-up  (a  detailed  treatment  of  this  function 
is  given  in  a  later  section) .  Since  the  run-up  occurs  in  the  close 
inshore  region  near  the  beach,  the  wave  angle  &  tends  toward  a 
small  angle,  so  that  dependence  of  the  vertical  force  and  pitch 
moment  in  terms  of  the  parameter  cos  y  can  be  neglected. 

11.  WIND  EFFECTS 

Forces  and  moments  due  to  wind  acting  on  the  above-water 
portions  of  the  hull  act  on  the  boat  in  the  horizontal  plane. 
Limited  data  on  wind  forces  exists  for  particular  ships,  but 
no  direct  information  is  available  for  the  present  assault  boat. 
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In  order  to  estimate  these  effects ,  it  is  necessary  to  determine 
the  relative  wind  speed  and  direction  with  respect  to  the  ship  in 
its  own  reference  frame.  Following  the  procedure  in  [26],  with 
the  wind  speed  Uw  and  direction  \^w  referred  to  the  inertial 

reference  frame  as  shown  in  Figure  9 ,  the  relative  wind  speed 
components  along  the  ship  axes  are 


URW  “  °w  °°S  ( V*1  '  U 
VRW  ’  UW  Sin  ‘V*1  '  V 


The  resultant  wind  is  then  given  by 

urw=  r 

at  the  angle  of  attack 
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relative  to  the  at,  as  shown  in  Figure  9, 
further  defined  uy 
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and  the  ship  yaw  orientation  are  measured  with  respect  to  the 
ine.  tial  reference  in  Figure  6  for  the  following  sea  condition. 
They  are  also  valid  when  measured  relative  to  the  same  reference 
for  determining  the  angle  i}<  for  the  case  where  the  boat  is 
heading  away  from  shore,  as  indicated  in  Figure  6,  due  to  the 
change  in  ship  axes. 

The  wind  force  in  the  direction  of  the  resultant  wind 
velocity  URW  is  given  in  [27]  as 


I)  -  _i  c  A  IT 

UW  2  u  proj  RW 


(163) 


where  p  is  air  density,  Cr  is  an  average  drag  coefficient,  and 
a 

\  .  is  the  silhouette  area  of  the  above-water  portion  of  the 

'  proi 

boat  projected  on  a  vertical  plane  normal  to  the  direction  of  the 
relative  wind  URVf.  The  quantity  Dw  in  Equation  (163)  is  computed, 
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and  assuming  a  small  lift  force  generated  by  the  boat,  the  compo¬ 
nents  of  Dw  are  resolved  along  the  boat  x-  and  y-axes.  With  A1 

the  effective  frontal  above-water  area  of  the  boat  hull  side,  and 
A2  for  the  lateral  above-water  area,  the  resulting  forces  are 

given  by 

p 

X (wind)  ~  ~  ~T  C d  (A1  icos  +  A2  lsin  I  >  urW  cos  1‘  (164) 

and 

p 

Y (wind)  =  ~T  CD  (A1  lcos  yl  +  A2  !sin  wl)  ujw  sin  u  (165) 

Since  the  boat  hull  is  almost  symmetrical  fore  and  aft, 
for  the  approximation  used  in  the  present  study,  the  above-water 
portion  may  also  be  viewed  in  that  way,  except  for  the  influence 
of  the  large  bow  ramp  and  the  pilot  house. 

The  moment  due  to  the  wind  depends  on  the  location  of  the 
(effective)  center  of  projected  area,  and  as  an  estimate  this  is 
represented  by 

N (wind)  =UI'ICD  <A1  lcos  "I  +  *2  !sin  » I  >  drw  ’ 

(1  -|  sin  p|)  sign  (sin  p)  (166) 

The  value  of  0.3  represents  the  limit  (as  a  fraction  of  ship 
length)  of  the  wind  pressure  centroid  as  the  wind  heading  changes, 
and  this  is  consistent  with  the  calculations  for  different 
vessels  presented  in  [25] . 

An  evaluation  of  the  numerical  parameters  appropriate  to 
the  LCM  (6)  hull  in  its  two  operating  modes  is  presented  in  [2], 

•  12.  BEACH  FORCES 

When  the  craft  lands  on  the  beach,  i.e.  reaches  a  point 
where  the  bow  portion  of  the  keel  touches  the  bottom,  it  is 
assumed  that  all  forward  motion  stops  and  the  beached  condition 
is  achieved  with  u  =  0.  In  that  case  forces  due  to  the  waves  in 
the  form  of  run-up  will  continue  to  act  on  the  craft,  especially 
in  the  lateral  plane,  where  the  yawing  moment  especially  requires 
attention  on  the  part  of  the  coxswain  to  provide  control  in  order 
to  maintain  the  craft  heading  normal  to  the  beach.  Due  to  the 
lack  of  knowledge  of  precise  soil  properties,  the  extent  of 
penetration  into  the  beach,  etc.  it  is  almost  impossible  to  deter¬ 
mine  the  forces  acting  on  the  craft  to  "hold"  it  on  the  beach  due 
to  the  interaction  between  the  beach  soil  material  and  the  craft. 
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most  significant  reactive  force  considered  is  that  in  the  axial 
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direction,  i.e.  the  force  that  must  be  overcome  by  the  ship 
propellers  in  order  to  "free"  it  from  the  beach.  When  this 
longitudinal  force  is  overcome  the  craft  can  move  out  into  deeper 
water  to  encounter  the  various  hydrodynamic  excitation  forces 
and  respond  to  them.  The  restraining  longitudinal  force  is 
assumed  to  be  equal  to  the  value  of  the  thrust  of  both  propellers, 
for  a  particular  value  of  propeller  rotational  speed  which  is 
selected  as  that  considered  representative  for  moving  a  landing 
craft  off  a  beach,  according  to  operational  experience  or  other 
guiding  dictates  of  the  simulator  training  officer.  Using  the 
expression  for  the  thrust  of  the  propellers  in  the  retraction 
mode,  the  beach  force  in  the  longitudinal  direction  is 

X (beach)  “  '  °-95  nt  (167> 

where  the  value  of  nfc  is  selected  for  a  particular  training 

mission.  This  same  force  expression  can  also  be  applied  when 
the  craft  is  encountering  the  beach  during  landing  in  order  to 
assist  in  bringing  it  to  a  final  stop. 

For  the  lateral  force  an  arbitrary  threshold  resisting 
value  can  be  selected  so  that  the  run-up  lateral  force  does  not 
cause  any  significant  sideward  motion.  The  yawing  action  must 
be  corrected  by  coxswain  control  via  differential  propeller  thrust, 
and  hence  no  yawing  moment  due  to  the  beach  is  considered.  Due 
to  small  roll  motions  anticipated  in  this  case  on  the  beach, 
just  as  in  the  case  of  run-up  effects  directly,  that  quantity  is 
also  neglected.  Since  the  vertical  force  and  pitching  moment  are 
due  to  hydrostatic  effects  of  the  run-up,  as  shown  in  Equations 
(156)  and  (157) ,  no  additional  beach  forces  on  these  modes  are 
considered,  since  the  resulting  motions  are  only  in  response  to 
the  run-up  flow,  which  occurs  in  a  time-varying  periodic  manner 
(see  discussion  in  a  later  section) . 
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B.  COMPUTATIONAL  CONSIDERATIONS 

The  equations  developed  in  the  preceding  analysis  are  to  be 
applied  to  a  simulator  that  represents  the  actual  motions  of  the 
craft  in  real  time.  Solutions  for  the  equations  must  be  obtained 
with  a  computer  that  would  also  be  able  to  represent  the  wave 
field  characteristics  in  proximity  to  the  boat,  as  well  as  the 
motions  of  the  craft  that  occur  simultaneously  due  to  the  wave- 
induced  excitation  and  the  coxswain-commanded  control  actions.  The 
equations  themselves  are  developed  in  a  mathematical  form,  but 
the  implementation  on  a  computer  is  a  separate  and  distinct  task. 
Examination  of  the  various  elements  entering  into  the  entire 
mathematical  model  of  the  craft  motion,  as  well  as  the  environ¬ 
mental  wave  field,  indicates  a  preference  for  a  particular  type 
of  computer  for  representing  each  element  in  the  overall  simula¬ 
tion. 

The  wave  generation  procedure  on  the  computer  is  best  carried 
out  by  means  of  a  digital  computer,  since  the  modifications 
necessary  to  account  for  the  influence  of  varying  depth  as  a 
function  of  time  as  the  craft  progresses  along  the  beach,  and 
the  transition  to  the  solitary  wave  form,  are  most  easily  carried 
out  using  that  type  of  computer.  The  logical  decisions;  the 
iterative  routine  to  compute  wave  number;  and  the  means  of 
computation  of  the  different  mathematical  functions  essential 
for  the  representation  of  a  wave  system  are  most  suitable  to  a 
digital  computer,  and  hence  that  is  the  recommended  procedure 
to  be  used.  Since  the  wave  excitation  forces  are  made  up  of 
functions  involving  wave  parameters,  as  well  as  similar  mathemati¬ 
cal  functions,  the  wave-induced  forces  and  moments  should  also  be 
computed  using  a  diaital  computer  technique.  The  computer 
operations  for  the  entire  simulation  assume  the  existence  of 
information  in  an  analog  domain,  which  is  converted  into  the 
digital  domain  by  use  of  A-D  converters,  so  that  the  digital 
computation  can  proceed.  The  output  of  the  diqital  program  is 
then  sent  through  D-A  converters  so  that  this  cata  is  then  in 
the  proper  form  for  direct  application  in  the  analog  representa¬ 
tions  of  the  dynamic  equations  of  motion, 

A  summary  of  the  final  equations  of  motion,  together  with 
functional  block  diagrams  and  associated  explanatory  discussion, 
is  presented  in  the  Appendix  of  this  report.  This  Appendix  also 
contains  the  flow  charts  for  computation  of  wave  characteristics , 
time  history  forms,  etc.  as  well  as  the  means  of  computing  the 
associated  wave  forces.  The  computational  technique  described 
above  is  directed  toward  use  of  both  analog  and  digital  computers, 
together  with  a  hybrid  interface.  The  use  of  a  hybrid  computer 
system  as  a  means  of  carrying  out  the  required  equation  solutions, 
etc.  is  not  proposed  as  the  only  possible  means  of  implementing 
the  present  mathematical  model.  However  it  is  the  type  of  system 
present  at  Qceanics,  and  it  can  provide  real  time  solutions  of  the 
dynamic  equations  of  motion  via  the  analog  portion  of  the  system. 
Assuming  that  a  general  purpose  digital  computer  with  proper 
capacity  were  available,  the  basic  technique  of  solution  indicated 
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in  the  flow  diagrams  could  also  be  implemented  on  that  type  of 
machine.  Questions  concerning  the  integration  method,  stability 
considerations,  sampling  times,  etc.  that  would  affect  the 
feasibility  of  carrying  out  such  computations  in  real  time  are 
considered  in  a  limited  fashion  in  a  later  portion  of  this  section. 

When  carrying  out  the  wave  field  computations,  the  concept 
proposed  is  to  be  able  to  evaluate  a  wave  elevation,  velocity,  etc. 
at  any  point  on  the  free  surface,  as  a  function  of  time,  independent 
of  the  location  of  the  craft.  Particular  locations  at  the  craft, 
as  it  progresses  on  toward  the  beach  or  as  it  retracts,  can  then 
be  considered  as  a  means  of  relating  the  environment  and  the  craft 
behavior.  This  point  of  view  also  allows  evaluation  of  the  various 
wave  forces,  by  referring  them  to  the  wave  as  measured  at  the  CG 
of  the  craft.  Similarly  it  allows  the  separate  computation  of 
oscillatory  waves,  solitary  waves,  and  run-up  elevation,  each  in 
the  appropriate  region  for  their  occurrence. 

The  spatial  and  time-varying  properties  of  the  oscillatory 
wave  and  the  corresponding  wave  forces  and  moments  are  properly 
simulated  as  shown  in  previous  sections  by  introducing  the  function 
f (t)  (defined  in  Equations  (123)  and  (132))  referenced  to  a  fixed 
point  offshore.  The  situation  is  different  with  respect  to 
solitary  wave  expressions  since  they,  by  their  hyperbolic  nature, 
represent  only  a  single  wave.  A  wave  train,  however,  can  be 
generated  in  a  simple  fashion  by  a  reinterpretation  of  the  argument 
g(t)  (defined  in  Equation  (140))  within  the  hyperbolic  functions. 

The  period  of  shoaling  breaker  waves  is  approximated  to  be  the 
period  of  maximum  energy  of  the  offshore  wave  spectrum,  T.  The 
crest  to  crest  distances  is  then  simply  cT  where  c  (the  wave 
propagation  speed)  varies  according  to  distance  offshore.  A 
train  of  solitary  waves  can  then  be  generated  if  the  argument 
g(t)  repeatedly  passes  through  zero  with  periodicity  cT.  The 
most  convenient  manner  of  implementing  this  in  a  digital  computer 
is  through  the  use  of  a  mod  function.  A  mod  function,  mod  (x,y) 
is  defined  as  x-(x/y]y  where  [x]  is  an  integer  with  magnitude  of 
not  more  than  x  and  sign  the  same  as  x.  It  is  easily  seen  that 
the  proper  function  is 

m ( t)  =  mod  | g  (t) ,cT  j  -  ~  (168) 

which  when  introduced  into  the  solitary  equations  for  wave  eleva¬ 
tion,  forces,  and  moments  will  give  the  correct  spatial  and  time- 
varying  properties. 

The  point  of  transition  between  oscillatory  waves  and  solitary 
waves  as  previously  discussed  is  known  once  the  offshore  wave 
spectrum  and  beach  slope  is  selected.  A  smooth  transition  is 
desirable  for  many  reasons,  but  particularly  with  respect  to 
forces  and  moments  so  that  unwarranted  transients  are  not  intro¬ 
duced  in  the  craft  motions.  A  smooth  transition  from  one  region 
to  the  other  can  be  most  easily  accomplished  if  a  transition 
zone  is  specified  in  which  the  wave  property  of  interest  from  the 
region  the  craft  is  leaving  is  decayed  out  while  the  similar 
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property  of  the  region  in  which  the  craft  is  entering  is  smoothly 
introduced.  The  width  of  the  transition  zone  is  chosen  to  be  about 
a  boat  length.  Thus  the  wave  elevation  in  the  transition  zone, 
upon  leaving  the  oscillatory  region,  could  be  represented  by 

(L“Sj  Sz 

n  (total)  =  — £ -  n(osc.)  +  -j-  n(sol.)  (169) 

where  s  is  the  distance  traversed  through  the  transition  zone, 
z 

where  the  transition  zone  extends  from  locations  of  length  j 

ahead  and  behind  the  tianci fion  point  (defined  by  relations  in 
Equation  (25)).  With  the  solitary  "-’ve  height  taken  equal  to  the 
significant  height  of  the  offshore  waves,  at  the  transition  point, 
the  height  grows  as  the  wave  propagates  into  shallower  water  in 
accordance  with  Greens's  law,  Equation  (27).  Breaking  occurs  when 
the  wave  height  and  the  local  water  depth  satisfy  Equation  (26) , 
and  this  point  can  be  established  by  initial  computations  when 
the  transition  point,  offshore  significant  height,  and  beach 
slope  are  given  or  established  (in  the  case  of  the  transition 
point,  in  terms  of  the  other  parameters  specified) . 

Another  problem  in  the  simulation  is  to  synchronize  the  run¬ 
up  with  the  solitary  waves  at  breaking.  Since  the  point  of 
breaking  is  known  beforehand,  through  the  use  of  the  mod  function 
the  times  at  which  the  solitary  wave  crests  arrive  at  the  breaker 
line  is  known.  These  are  the  times  at  which  the  run-up  function 
is  initiated. 

The  run-up  flow  field  is  established  following  breaking,  and 
it  is  assumed  to  have  a  time  extent  of  ,  with  T  the  period  of 

maximum  spectral  energy  offshore.  This  is  followed  by  a  time 
T 

interval  of  ^  in  which  no  elevation  occurs,  and  then  the  system 

repeats  over  and  over  again  as  a  time  history,  similar  to  a  series 
of  square  waves.  This  idealized  model  is  assumed  to  have  a 
propagation  speed  cR  defined  by  Equation  (32) ,  and  the  spatial 

and  temporal  variation  is  expressed  in  terms  of  the  function 

s  =  xR  -  cRt  (170) 

where  xD  is  the  longitudinal  distance  measured  from  the  breaking 

point,  along  the  x  -axis,  and  the  time  t  in  Equation  (170)  is 

measured  from  the  instant  of  wave  breaking  in  each  case.  To 
relate  this  result  to  other  parameters  developed  ii  the  analysis, 
the  relation 

XR  ‘  (VV  >  +  5  (171) 

is  defined,  where  x  is  the  inertial  frame  coordinate  (along 

Ou 
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xQ-axis)  at  which  breaking  occurs,  xQ  is  the  reference  position  of 

the  ship  CG  (or  any  space  points  at  which  the  run-up  form  is 
desired  to  be  determined) ,  and  5  is  the  "dummy"  variable  along  the 
ship  hull  (measured  from  the  origin  at  the  CG)  that  allows  deter¬ 
mination  of  the  relative  position  of  the  run-up  elevation.  The 
run-up  is  then  functionally  represented  by  the  expression 


R(xo,t)  =  R 


1 


(s,cRT) 


=  R  fR(«,t)  ,  (172) 

the  unit  step  function,  where  the  function  fRU,t) 

in  Equations  (156)  and  (157)  is  then  defined. 

Another  remaining  consideration  is  the  question  of  sampling 
rate  in  real  time  at  which  the  wave  forces  and  moments  have  to  be 
evaluated  and  inputted  to  the  dynamic  system  implemented  on  the 
analog  portion  of  the  computers.  The  sampling  theorem  requires 
(see  e.g.  [28])  the  sampling  rate  to  be  at  least  twice  the 
highest  frequency  of  interest  present  in  the  dynamics  to  avoid 
aliasing.  Since  the  bandwidth  is  about  2.5  times  the  natural 
frequency,  for  this  craft  the  sampling  rate  should  be  at  least 
5  times  the  highest  natural  frequency.  This  sampling  rate  should 
be  considered  to  be  only  a  lower  bound  and  a  higher  sampling  rate, 
perhaps  10  or  more  times  the  highest  natural  frequency  is  more 
desirable. 


with  1 


If  the  overall  simulation  is  to  be  carried  out  in  real  time 
by  means  of  a  digital  computer,  it  is  necessary  to  choose  an  appro¬ 
priate  integration  scheme  for  that  purpose,  with  the  various  require¬ 
ments  of  the  wave  field,  forces,  etc.  computations  coordinated 
consistently  within  the  solution  technique.  A  required  criterion 
of  a  suitable  integration  technique  is  efficiency,  i.e.  that  it 
accomplish  the  required  integration  in  a  minimum  of  time  so  that 
real  time  simulation  be  accomplished.  Another  requirement  is 
accuracy,  but  accuracy  is  often  indirectly  related  to  efficiency 
since  a  reduction  in  the  integration  interval  often  produces  a 
concurrent  increase  in  accuracy.  A  third  criterion  is  stability, 
where  the  integration  method  does  not  permit  any  accumulated 
errors  to  grow  or  oscillate  without  bound  over  the  length  of 
time  of  particular  simulation  applications.  All  of  these  factors 
are  related,  in  very  intricate  ways,  and  there  is  no  single 
guidance  to  indicate  a  choice  of  any  single  method  for  the  present 
problem  (or  for  that  matter,  for  any  arbitrary  real-time  vehicle 
simulation).  However,  while  there  are  many  integration  schemes 
to  choose  from,  .'or  a  general  systems  simulation  where  speed  and 
accuracy  is  of  paramount  importance  and  numerical  stability  is 
assured,  the  Runge-Kutta-Merson  method  [29]  is  recommended  on 
the  basis  of  a  comparative  study  of  the  performance  of  this  method 
with  the  other  common  methods  [30] . 
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SECTION  IV 
CONCLUSIONS 

The  six  degree  of  freedom  mathematical  model  for  an  assault 
boat  developed  herein  is  concluded,  on  the  basis  of  the  experi¬ 
ments  presented  in  [2] ,  to  be  an  accurate  representation  of  the 
boat  motions  from  both  the  hydrostatic  and  hydrodynamic  effects 
of  body-fluid  interaction  in  the  presence  of  waves.  The  waves 
are  considered  to  have  been  generated  by  an  offshore  random  sea 
and  are  progressing  towards  the  shoreline.  Account  is  taken  of 
the  change  in  wave  character  and  direction  resulting  from  the 
influence  of  a  sloping  beach  to  the  time  of  breaking,  which  is 
followed  by  run-up  onto  the  beach.  Effects  of  run-up  and  other 
causes  of  boat  motion,  namely,  propulsion  and  wind  effects  and 
the  coxswain  control  commands  are  detailed. 

It  is  concluded  that  this  model,  together  with  the  computa¬ 
tional  techniques  discussed  herein,  serves  as  an  adequate  basis 
for  the  development  of  a  real  time  simulator  for  the  purposes 
of  coxswain  training. 
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Figure  2.  Variation  of  wave  characteristics,  relative  to  deep  water 

values,  as  function  of  ratio  of  depth  to  deep  water  wavelength 
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Figure  3.  Wave  refraction  diagram  for  waves  approaching 
a  uniformly  sloping  beach,  8Q  *  45° 
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"offshore " 


Figure  4.  Waves  approaching  beach  over  varying  bottom  contour, 
with  spatial  reference  point  locations 
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Figure  5.  Sea  State  4  Wave  Spectrum  Decomposition 
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shoreline 


Figure  6.  Relation  of  ship  axes,  shore  position,  and 
inertial  reference  frame 
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APPENDIX  A 


EQUATION  SUMMARY  AND  SOLUTION  PROCEDURES 


The  final  equations  of  motion  are  represented  below,  in 
general  form  and/or  with  reference  to  other  expressions  within 
this  report,  for  operating  conditions  corresponding  to  approach¬ 
ing  (heavy  configuration)  or  retracting  (light  configuration) 
from  a  sloping  beach  through  waves.  The  neglect  of  coupling 
between  planes  of  motion,  the  reference  to  axes  appropriate  to 
seakeeping  analyses,  and  other  simplifications  reduce  or  delete 
various  nonlineiu:  effects,  thereby  allowing  ease  of  computational 
requirements  for  solution  of  the  equations.  The  compiled  equa¬ 
tions  are  presented  for  each  degree  of  freedom,  with  a  separate 
discussion  following  the  presentation  to  indicate  the  various 
procedures  (in  flow  chart  form)  for  computing  the  wave  properties, 
forces,  etc.  together  with  their  use  in  the  equations.  The 
procedures  for  solution  of  the  equations  are  also  outlined  in  a 
flow  diagram,  and  reference  is  given  to  particular  equations  in 
the  main  text  where  various  parameters,  operations,  etc.  are 
defined  and/or  derived.  Amplifying  discussion  is  given,  where 
required,  to  explain  some  of  the  procedures  in  the  computations. 
The  numerical  values  appropriate  to  various  parameters  in  the 
equations  are  given  in  [2]  or  are  pre-selected  for  each  particu¬ 
lar  problem  (such  as  sea  conditions,  initial  conditions,  etc.). 


Accelerations 


u 


ra-X. 


u 


<  X(hull)  +  X (props)  +  X (rudders)  +  X(wind) 


+  X (waves)  +X (beach)  +  ravr 


(A-l) 


v  “  m-Y»  Vluv  +  <Y2~w)ur  +  Y  (props)  *  Y  (rudders) 


+  Y(wind)  +  Y (waves) 


(A-2) 


z  = 


m+2£A 


33 


-2£A'33U0 


2£N*  (z+U0) 

z  z 


2pg*Bz  +  Z (waves) 


(A- 3) 
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♦  = 


Jxx”Kp 


K  i-  W|GMU+  K(waveJ 


(A-4) 


J 


V  f  i3A33 


'e'  =  7~ . f  T7"  lAA33{u*+u2e)  -  I  -  !  P9B^e 


r  = 


+  M . 
1 


l  waves ; 
(  . 


X  -N* 
zz  r 


fA-5^ 


<  Nlur  +  N2rlrf  +  N3uv  +  N (props)  +  N  (rudders) 


+  N  (control)  +  N(wind)  +  N (waves) 
Velocities  and  Displacements 


(A-6) 


Velocities  and  displacements  are  determined  by  integration, 
leading  to  typical  results  such  as 

ft 


and 


u  = 


*  * 


u  dt  ,  etc. 


r  dt  +  \i>  (0) 


(A-7) 


(A-8) 


Kinematics 

In  addition  to  the  above  equations,  kinematic  relations 
providing  craft  position  in  space  with  respect  to  the  inertial 
reference  frame  are  necessary.  These  relations  are  obtained 
from  Equations  (52)  and  (53)  and  are  expressed  as 


xo  = 


ft 


u  dt  *  cos  \f)  - 


v  dt  •  sin  ip 


(A-9) 


r  t 


v  dt  •  cos  ip  + 


J  o 


u  dt  •  sin  ^ 


(A- 10) 


for  the  condition  when  proceeding  toward  shore  (heavy  configura 
tion) ,  and 
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(A-ll) 

(A-12) 

when  going  seaward  from  the  beach  (light  configuration) . 

Hull  Forces 

X(hull)  is  9iven  by 


x_  = 


ft 

v  dt  •  sin  ip  - 
Jo  J 


u  dt  •  cos  \p 


yn  -  - 


v  dt  •  cos  \p  - 


u  dt  •  sin  ip 


X(huU>  =  Xlulvl  +  X2vr  <A-13) 

for  the  heavy  configuration,  and  by 

X<huUi  -  X2vr  +  X3vS  (A-14) 

for  the  light  configuration.  All  the  hull  forces  for  other 
modes  of  motion  are  expressed  in  the  acceleration  expressions  in 
Equations  (A~l)  -  (A-6) . 

Propeller  Forces 

X (props)  is  rePresented  by  Equations  (88) 

and  (90)  for  the  heavy  configuration,  and 
by  Equations  (91)  and  (92)  for  the  light 
configuration.  (A~15) 

The  rotational  speeds  of  the  port  and  starboard  propellers 

are  represented  by  n  and  n  (in  rps)  and  the  operating  propeller 

P  3 

speed  used  in  the  propeller  force  equations  is  defined  as 

n  =  ~  (n  +n  )  (A-16) 

O  2  p  S 

Y (props)  rePresented  by  Equation  (84) 

for  the  heavy  configuration,  and  by  Equation 

(86)  for  the  light  configuration.  (A-17) 


N (props)  5^ven  by  Equation  (85)  for  the 

heavy  configuration,  and  by  Equation  (87) 

for  the  light  configuration.  (A-18) 
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Rudder  Forces 

X{ rudders)  is  glven  by 

"(rudders)  =  '  a-8u2s2 

x (rudders)  is  9iven  b* 

Y (rudders)  *  2'71  u3{  (heavy> 

"(rudders)  "  °'47  “2s  ai*ht) 

"(rudders)  is  ®iven  b" 

"(rudders)  =  '  56-63  “3{  (heavx) 

"(rudders)  '  12'28  “3{  <llc-'bt> 

where  s  is  expressed  in  terms  of  radians. 

Control  Moment 

The  yawing  control  moment  due  to  differential 
propeller  thrust  is  given  by  Equation  (110) , 
using  the  propeller  thrust  representation  for 
either  the  heavy  or  light  configuration  given 
in  Equations  (104)  and  (105) ,  as  the  particu¬ 
lar  case  requires. 


Wind  Forces 


X(wind)  is  given  by  Eciuation  (164);  Y  (wind) 
is  given  by  Equation  (165);  and  N(win<j)  is 

given  by  Equation  (166).  These  terms  are 
expressed  in  terms  of  parameters  defined  in 
Equations  (158)  -  (161) . 


Beach  Forces 


X  (beach)  9;‘-ven  ky  -Equation  (167)  and  the 

associated  discussion  related  to  that  equa¬ 
tion  in  the  text  of  the  report. 
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Wave  Forces 


X  (waves)  rePreseTVfce<^  hY  Equation  (127) , 

followed  by  a  transition  region  (involving 
"blending’*  of  force  representations,  to  be 
discussed  when  considering  wave  effects)  and 
then  the  representation  given  by  Equation 
(142).  This  is  followed  by  Equation  (152) 
when  the  craft  is  in  the  run-up  region 
following  breaking  (a  discussion  of  the 
particular  region  and  the  timing  operations 
to  apply  to  the  run-up  are  given  in  this 
Appendix) .  Another  part  of  the  longitudinal 
wave  force  is  given  by  Equation  (134) ,  which 
is  added  to  the  above  when  the  craft  is  in 
the  region  of  oscillatory-type  waves,  in 
transition,  and  when  in  the  solitary  wave 
region,  where  all  of  these  expressions  are 
valid  for  the  heavy  configuration  approaching 
the  shoreline.  For  the  light  configuration 
the  wave  force  is  represented  by  Equation 
(153)  in  the  run-up  region*  followed  by 
Equation  (142)  with  the  sign  of  u  changed 
and  then  through  transition  into  Equation 
(127)  with  the  sign  of  u  changed.  To  this 
is  added  Equation  (135)  in  the  region  of 
solitary  waves  through  transition  into  the 
oscillatory-type  wave  region. 


Y (waves)  9*ven  Equation  (128),  followed 

by  a  transition  into  the  expression  in  Equa¬ 
tion  (148) .  This  is  then  followed  by  Equa¬ 
tion  (154)  in  the  run-up  region,  where  all  of 
those  expressions  are  valid  for  the  heavy 
configuration.  For  the  light  configuration 
the  lateral  wave  force  is  given  by  Equation 
(154) ,  followed  by  Equation  (147) ,  and  then 
Equation  (128) ,  where  the.  sign  of  u  is 
changed  throughout  these  expressions  to 
reflect  the  altered  orientation  of  the  craft. 


'J  (waves)  9^ven  bY  the  sum  °f  Equations 

(124)  and  (125) ,  followed  by  a  transition 
to  Equations  (144)  and  (146) .  Then  following 
breaking  by  Equation  (156) ,  with  the  defini¬ 
tions  in  Equations  (170)  -  (172) ,  where  t'e 
results  are  for  the  heavy  configuration.  For 
the  light  configuration  the  v<-'  ticc.l  wave 
force  is  given  by  Equation  u.5  ,  followed 

by  the  sum  of  Equations  (144)  and  (146) ,  and 


(A-27) 


(A-28) 
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then  by  the  euro  of  Equations  (124)  and 
(125) ,  where  the  signs  of  the  terms  in 
these  last  equations  (Equations  (144) , 
(146),  (124)  and  (125))  are  changed. 

In  addition  the  sign  of  u  is  also 
changed  throughout  these  expressions. 


K (waves)  “or  heavy  configuration  is 

given  by  Equation  (131) ,  followed  by 
Equation  (150)  after  transition  into  the 
solitary  wave  region.  After  breaking  no 
other  roll  moment  is  assumed  to  act  on 
the  craft  due  to  any  wave  action.  For  the 
light  configuration  K(waves)  is  represented 

in  the  solitary  wave  region  by  Equation 
(150),  and  then  by  Equation  (131),  where 
the  sign  of  u  is  changed  in  these  express¬ 
ions  for  this  case. 


M (waves)  ^or  heavy  configuration  is 

given  by  Equation  (126) ,  followed  by 
transition  into  the  sum  of  Equations 
(145)  and  (147).  This  is  then  followed 
by  Equation  (157)  in  the  run-up  region. 
For  the  light  configuration  M(waves) 

is  given  by  Equation  (157) ,  followed 
by  the  sum  of  Equations  (145)  and  (147) 
in  the  solitary  wave  region,  and  then 
by  Equation  (126)  after  the  transition 
zone.  The  sign  of  u  is  changed  in  all 
the  expressions  above  for  the  light 
configuration. 


N (waves)  for  heav*  configuration  is 

given  by  Equation  (130),  followed  by 
transition  into  Equation  (149) ,  and  then 
following  breaking  by  Equation  (155)  in  the 
run-up  region.  For  the  light  configuration 
N (waves)  9-*-ven  by  Equation  (155), 

followed  by  Equations  (149)  and  then 
Equation  (130)  after  transition.  The 
signs  of  the  terms  in  Equations  (149)  and 
(130)  are  changed,  and  also  the  sign  of  u 
wherever  it  appears  in  these  equations 
for  the  light  configuration. 
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A  flow  diagram  for  solution  of  the  equations  developed  in 
this  study  is  given  in  Figure  A-l.  All  of  the  constituent 
elements  presented  in  this  Appendix  are  shown  there,  together 
with  reference  to  particular  equations,  descriptive  paragraphs, 
or  other  explanatory  sources  in  the  Appendix.  This  will  relate 
all  of  the  terms  and  operations  required  in  carrying  out  a 
computer  simulation  using  these  equations. 

When  treating  the  waves,  and  from  them  the  wave  forces,  a 
specific  sequence  of  operations  is  required  to  represent  the  wave 
form  and  forces.  An  illustration  of  the  procedures,  with  ref¬ 
erence  to  specific  equations  in  the  body  of  the  report,  is  given 
in  Figure  A-2  for  the  oscillatory-type  waves  and  in  Figure  A-3 
for  solitary  waves.  Since  the  wave  forces  and  moments  contain 
the  same  functional  expressions  as  the  wave  form  expressions, 
their  computation  readily  follows  after  the  waves  are  determined. 

The  separate  regions  where  different  types  of  waves  are 
present  are  defined  by  the  following: 

Oscillatory  region,  h  >  .04  Xq  +  1.35  Hq  (A-33) 

Solitary  region,  h  <  .04  XQ  +  1.35  Hq  (A-34) 

H  5 

Run-up  region,  when  g-  >  j  (A- 35) 

where  h  is  local  depth,  Hq  is  the  offshore  significant  wave 

height,  Aq  is  the  representative  offshore  wavelength  (taken  from 

period  of  maximum  spectral  energy) ,  and  H  is  the  solitary  wave 
height  (computed  continuously,  although  not  used  beyond  breaking). 
Methods  of  determining  the  breaking  location  as  well  as  the  timing 
sequence  for  initiating  run-up  functions  are  discussed  in  a  later 
section  of  this  Appendix. 

The  method  of  selecting  the  appropriate  wave  forces  for 
inclusion  in  the  flow  diagram  of  Figure  A-l  is  illustrated  in 
Figure  A-4.  Within  that  diagram  the  operation  of  "sum  all  i’s" 
means  the  summation  of  all  the  separate  contributions  of  the  ele¬ 
mental  waves  that  make  up  the  offshore  spectrum  (i  =  1 *  8  chosen 
for  the  present  case) ,  and  similarly  for  the  wave  forces  and  moments 
from  each  of  the  component  terms.  The  "blend"  operation  in  that 
figure  refers  to  the  transition  treatment  for  the  various  wave 
forces  as  well  as  the  wave  form  itself,  when  combining  oscillatory 
and  solitary  wave  terms  in  the  transition  region.  This  region 
is  assumed  to  be  of  length  L,  the  boat  length,  and  extends 

±  j  in  either  side  of  the  transition  loc;  tion  given  by  Equation 

(25)  in  the  text  (or  by  the  equality  of  the  relations  in 
Equations  (A-33)  and  (A-34) .  For  -he  case  of  the  wave  motion  as 
an  example,  when  approaching  -he  beach  seaward,  the  blended 

wave  elevation  is  represented  by 
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n 


,  x 

nosc.  L  nsol. 


(A- 36) 


where  n  can  be  held  constant  at  the  last  calculated  value 
osc. 

before  entering  the  transition  region.  The  quantity  x  \s  defined 
as 

(A-37) 

which  is  the  distance  penetrated  into  the  transition  zone,  where 
xfc  is  the  distance  along  the  xQ-axis  from  the  origin  to  the 

transition  location  (Equation  (25)).  When  coming  off  the  beach 
and  proceeding  seaward,  the  last  solitary  wave  values  are  held 
constant  in  the  computation  of  Equation  (A-36) ,  and  the  value 
of  x  is  given  by 


(A-38) 


The  location  where  solitary  waves  break  and  run-up  onto  the 
beach  can  be  calculated  directly  from  the  selected  offshore 
wave  spectrum  characteristics  and  beach  slope.  It  has  been 
previously  shown.  See  Equation  (25)  that  these  same  parameters 
also  determine  the  transition  point  at  which  the  oscillatory 
waves  become  solitary  waves.  Denoting  the  depth  at  which 
transition  occurs  as  hfc,  and  assuming  the  initial  solitary  wave 

height  at  transition  to  be  equal  to  the  significant  wave  height 
of  the  offshore  wave  spectrum,  Hq,  the  height  of  a  solitary  wave 


as  it  progresses  toward  the  beach  is  governed  by  Equation  (27) 
with  =  ht  and  =  Hq  until  a  wave  height  to  local  depth 

ratio  of  .714  (=  y)  is  obtained  (Equation  (26)).  Combining  these 


equations  and  solving  for  the  depth  at  which  breaking  occurs 


which,  in  turn,  knowing  the  beach  slope,  determines  the  distance 
Xq  .  It  is  important  to  realize  that  breaking  of  a  soJitary 
b 

wave  may  not  occur  at  the  time  the  landing  craft  reaches  the 
breaking  point  and,  in  fact,  the  coxswain  should  attempt  to  avoid 
this  situation  by  adjusting  the  boat  speed  so  that  the  boat 
follows  behind  a  solitary  wave  and  its  resulting  run-up  onto 
the  beach. 


The  run-up  is  synchronized  in  time  with  the  arrival  of  a 


solitary  wave  at  the  breaking  point, 


,  by  simply  replacing 
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T 

time,  t/  in  all  run-up  functions  with  t+  j  .  This  is  due  to  the 

T  ^ 

offset  of  j  seconds  introduced  in  Equation  (168)  ho  achieve  the 

proper  time  characteristics  of  the  hyperbolic  functions. 
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Figure  A-l.  Flow 
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?r  m  for  solution  cf  eruptions  cf  motion 
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